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Abstract 

We construct a family of BC„-symmetric biorthogonal abelian functions generalizing Koornwinder's 
orthogonal polynomials, and prove a number of their properties, most notably analogues of Macdonald's 
conjectures. The construction is based on a direct construction for a special case generalizing Okounkov's 
interpolation polynomials. We show that these interpolation functions satisfy a collection of generalized 
hypergeometric identities, including new multivariate elliptic analogues of Jackson's summation and Bailey's 
transformation. 
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1 Introduction 



In |16j . we considered two families of _BC n -symmetric Laurent polynomials: the multivariate orthogonal polyno- 
mials introduced by Koornwinder [Sj, and the "interpolation Macdonald polynomials" of Okounkov related 
by Okounkov's "binomial formula" . One of the results of that paper was a multivariate analogue of Jackson's 
summation, which had originally been conjectured by Warnaar |29|. In fact, that summation was a limiting 
case of Warnaar's conjecture, an identity of elliptic hypergeometric series. This suggested that the theory of 
[T?)j should be extensible to the elliptic case; this extension is the topic of the present work. 

The main difficulty in constructing this extension is the result of Okounkov ^3] which roughly speaking said 
that the interpolation Macdonald polynomials could not be generalized; more precisely, he defined a general 
notion of interpolation polynomial, and showed that the Macdonald-type polynomials form a component of the 
corresponding moduli space. Thus in order to extend our results to the elliptic level, it would be necessary to 
generalize this notion of interpolation polynomial. 

Another indication that an elliptic extension should exist was the existence of elliptic analogues [57||2S|, both 
discrete and continuous, of the Askey- Wilson orthogonal polynomials [2] • The difference from the Askey- Wilson 
theory (already seen at the trigonometric level in work of Rahman |14j for continuous biorthogonality, and later 
work by Wilson [3U| for discrete biorthogonality) was two-fold: first, rather than being orthogonal, the functions 
were merely biorthogonal , and second, the constraint of being monic of specified degree became a constraint 
on allowable poles. Extending this prescription directly to the multivariate level ran into the difficulty that 
divisors in multiple variables are much more complicated than in the univariate case. However, if we first clear 
the poles, we find that the constraint on poles becomes a requirement that the function vanish at appropriate 
points. 

We thus arrive at the notion of "balanced" interpolation polynomial; this differs from Okounkov's notion 
in that, rather than have one collection of vanishing conditions together with an assumption of triangularity 
on monomials, balanced interpolation polynomials satisfy two complementary collections of vanishing condi- 
tions. (This definition was independently discovered by Coskun and Gustafson 4 , who also obtained many 
of the results of Section in particular Theorems 14.11 and 14.91 below.) This allowed the affine symmetry of 
ordinary interpolation polynomials to be extended to a projective symmetry, and as a result allowed the desired 
generalization to elliptic interpolation polynomials. Moreover, two new special cases arise: first, by matching 
up the vanishing conditions, one obtains a collection of delta functions; more subtly, another special case ap- 
pears in which the interpolation polynomials factor completely. These special cases then play the role that the 
Macdonald polynomials did in , namely that of instances that can be solved directly, thus giving boundary 
conditions for recurrences. 

In particular, we have been able to generalize many of the results of |16j , especially concerning the Koorn- 
winder polynomials, which we generalize to a family of biorthogonal abelian functions. A closely related family 
was studied in using a certain contour integral identity and a related integral operator; in particular, it was 
shown there that the biorthogonal functions satisfy an analogue of Macdonald's normalization and evaluation 
conjectures. The analogue of the remaining conjecture (referred to here as "evaluation symmetry" ) requires a 
more in depth understanding of interpolation polynomials, and will be proved in Theorem l5.4l below. In contrast 
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to |15j . our present approach is essentially algebraic in nature; although to begin with we use the theory of 
theta functions over C, we will eventually see that all of our results can be given purely algebraic proofs. 

We begin in Section 2 by giving the definition and basic properties of balanced interpolation polynomials, 
in particular considering four of the main cases satisfying extra vanishing properties (the "perfect" case); the 
fifth ("elliptic") case is the subject of the remainder of the paper. 

Section 3 begins our treatment of the elliptic case by defining a family of £>C„-symmetric theta functions, 
parametrized by partitions contained in a rectangle, defined via vanishing conditions. As in |16|. the key 
property of these functions is a difference equation, from which the extra vanishing property follows. We also 
define a corresponding family of abclian functions, this time indexed by all partitions of length at most n. 

Evaluating an interpolation abelian function at a partition gives a generalized binomial coefficient; we study 
these coefficients in Section 4. In particular, we obtain a pair of identities generalizing Jackson's summation 
and Bailey's transformation respectively (sec 6 for the univariate elliptic case); each identity involves a sum 
over partitions contained in a skew Young diagram of a product of binomial coefficients. Using these, we derive 
some unexpected symmetries of binomial coefficients, relating a coefficient to coefficients with conjugated or 
complemented partitions. This then enables us to prove a number of properties of interpolation functions, 
including a connection coefficient formula, a branching rule, a Pieri identity, and a Cauchy identity. 

Section 5 introduces the biorthogonal functions, defined via an expansion in interpolation functions, and 
shows how they relate to the functions considered in Using the interpolation function identities, we prove 
a number of identities for the biorthogonal functions; in addition to the analogue of Macdonald's conjectures 
for Koornwinder polynomials, we obtain connection coefficients, a quasi-branching rule, a quasi-Pieri identity, 
and a Cauchy identity. 

In Section 6, we depart from the analytic approach via theta functions, and show that, in fact, the inter- 
polation and biorthogonal functions can be defined purely algebraically (and thus, over the complex numbers, 
are invariant under transformations of the modular parameter) . We sketch a purely algebraic derivation of the 
identities of Sections 4 and 5; the key step being to construct the algebraic analogue of the difference operator 
of Section 3. 

In Section 7, we use our results on interpolation functions to define a family of perfect bigrids based on 
elliptic curves. Since this family is an open subset of the space of perfect bigrids, we in particular find that 
we cannot add any parameters to our elliptic theory, without further extending the notion of interpolation 
polynomial. We then consider how the bigrids in this family can degenerate, obtaining four general classes of 
degeneration (corresponding to the Kodaira symbols I\, I2, II, and III; these in turn split somewhat further 
over a non-algebraically-closed field), including ordinary hypergeometric series (// and some instances of III) 
and basic hypergeometric series {I\ and some instances of I2) as well as some apparently unstudied cases. 

Section 8 shows how the present results relate to the results of ^B] . The identities of that paper involved 
three different functions of skew Young diagrams: binomial coefficients, "inverse" binomial coefficients, and 
principal specializations of skew Macdonald polynomials. We show that each of these can be obtained as a limit 
of elliptic binomial coefficients, and thus in particular prove that the Koornwinder polynomials are a limiting 
case of the biorthogonal functions. 

Finally, in Section 9, we consider some open problems relating to interpolation and biorthogonal functions. 
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Notation 

As in [lfij . we use the notations of ^Ij f° r partitions, with the following additions. If A C m", then to" — X 
denotes the "complementary" partition defined by 

{m n ~\)i =m- X n+ i-i, l<i<n (1.1) 

In addition, if ^(A) < n, we define to" + A to be the partition 

(m n + A)i = m + Xi, l<i<n. (1.2) 

Similarly, if Ai < to, we define to" • A to be the partition such that 




(m n -X)i = i (1.3) 



(Note that this is specifically not defined for Ai > to; it is a concatenation operation, not to be confused with 
a union operator.) We also define a relation -< m on partitions such that K -< m X if and only if 

kCActo"+k (1.4) 

for sufficiently large n; similarly, k -<' m X iff k! -< m A'. (If m = 1 in either case, the subscript will be omitted.) 

We will also need the following "elliptic" analogues of g-symbols and generalizations. Let p be a complex 
number such that < \p\ < 1. Then we define 

9(x;p):=l[(l-p k x)(l-p k + 1 /x) (1.5) 

0<fc 

e(x;q;p) m := ]J 0(q'x;p) (1.6) 

0<j <m 

In each case (and similarly for the C symbols (but not the A symbols) below), we follow the standard convention 
that the presence of multiple arguments before the (first) semicolon indicates a product; thus for instance 
9(ax ,p) — 6(ax;p)9(a/x;p). We also need elliptic analogues of the C symbols of |16) : 

C° x (x;q,t;p) := [] fftf-H^xip) (1.7) 

(i,j)6A 

C;(x;q,t;p):= JJ 6( q x *-H x '-*x;p) (1.8) 

C+(x;q,t;p):= [] B^-H^-^p). (1.9) 

Note that in the limit p = we recover the symbols of JHI ; similarly, the analogous symbols of ^3 can each be 
expressed as products of two of our present symbols (with q and p switched in one symbol). The transformations 
of |16| carry over to the present case, via the same arguments. 

Two combinations of C symbols are of particular importance, and are thus given their own notations. For 
any nonzero complex number a and any finite sequence of nonzero complex numbers . . . 6,; . . . , we define two A 
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symbols as follows. 

A x (a\ . . . h . . . ; q, t;p) := A A (a| ...b t ...;q, t;p) — ^ + — . (1.11) 

C\ (pq,t; q,t;p)C\ {a,pqa/t; q,t;p) 

(Here, as in [16|. 2A 2 denotes the partition such that (2A 2 ).; = 2(A^/ 2 ])-) Note that although A essentially 
follows the standard multiple argument convention, A very much does not. We note the following transforma- 
tions: 

A x (a/qt\...b l ...;q,t;p) (1.12) 

j.2n— 2 j-n—lfo 

A A(^-| • ■ ■ ■ ■ -,t U , q~ m ,pqt n -\pq/q m t; q,t;p) (1.13) 



A v (o| 


...6, 


...;l/t,l/q;p) 


A m „_ 


-A(a|. 


..bi...;q,t;p) 


A ro r, 


,(a|.. 


.bi...;q,t;p) 


A m ti_| 


-a(o| ■ 


..bi...;q,t;p) 


A m r. 


,(a|.. 


■ k.. -;q,t;p) 


A m „ 


•a(o|. 


..b t ...;q,t;p) 


A w 


a(a|.. 


■ bi...;q,t;p) 



q 2m a q"°a 

A x (q 2m a\ . . . q m bi t n ,pqt n -\q m a/t n -\pqq m a/t n - q, t;p) (1.14) 

A A (a/i 2 "| . . . bi/1* q- m lP q/q m t, q m a/t n '\pqq m a/t n - q, t;p), (1.15) 

with corresponding transformations for A ; we also note that 

A m r,(a| ...bi.. .;q,t;p) = lim A^„(a| ...bi.. . ,q m pqa/x,pqa/xt n ,pqax,pqt n ~ 1 x/q m ;q,t;p) (1.16) 

A meromorphic (p-)theta junction is a meromorphic function / on C* (the multiplicative group of C; i.e., / is 
meromorphic away from and oo) such that 

f(x) = ax n f(px) (1.17) 

for suitable a <E C*, n 6 Z; without the adjective "meromorphic", theta functions are assumed holomorphic. 
Thus in particular the function 8(x;p) is a theta function, satisfying 

0( X ;p) = -x6(l/x;p) = -x9(px;p). (1.18) 

An elliptic function is a meromorphic function on C* invariant under multiplication by p; in other words, a 
meromorphic function on the quotient group C* / (p) (which is a compact 1-dimensional complex Lie group, or 
in other words an elliptic curve). 

We also have the following classes of multivariate functions. 

Definition 1. A meromorphic BC n - symmetric theta function of degree m is a meromorphic function / on 
(C*) n such that 

f(xi, . . . , x n ) is invariant under permutations of its arguments. 
f(x\, . . . , x n ) is invariant under Xi i— > 1/xi for each i. 
f(px 1 ,x 2 ,..-,x n ) = (l/pxj) m f(x 1 ,x 2 , ■ ■ ■ ,x n ). 
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A BC' n - symmetric (p-)abelian function is a meromorphic _BC„-symmctric thcta function of degree 0. 

Remark 1. Again, i3C n -symmetric theta functions are assumed holomorphic unless specifically labeled "mero- 
morphic" . 

Remark 2. It is not entirely clear to which root system these functions are truly attached. Certainly, the 
biorthogonal functions considered below generalize the Koornwinder polynomials, traditionally associated to 
the root system BC n , but there is also good reason to associate this class of theta functions to the affine root 

( 21 it^i _____ 

system of type A n (see, for instance, the work of Looijenga an d Saito on the invariant theory of abclian 
varieties associated to (extended) affine root systems.) In the absence of a general theory, we have chosen the 
more familiar label. 

The canonical example of this is the function Y[i<i<n^( ux i'yP) f° r u £ C*; indeed, the graded algebra of 
-BC„-symmetric theta functions is generated by such functions. 

We define a difference operator on (meromorphic) _BC n -symmetric theta functions that will play a crucial 
role in the sequel; this is the analogue of the difference operator of JHI i an d also appeared prominently in JS] . 

Definition 2. Let n be a nonnegative integer, and let a, 6, c, d G C* be arbitrary parameters. Define a 
difference operator D^ n ' (a,b, c, d; q,t;p) acting on i3C„-symmetric meromorphic functions on (C*) n by: 

(£> (n) (a,6, c, d;q,t;p)f)(x 1 ,...x n ) 

a£{±l}» KKn U[ - X i > P > l<i<j<n 1 * X J 1 P) 

Proposition 1.1. For each nonnegative integer m, if q m t n ~~ 1 abcd — p, then D^ n '(a,b,c,d:,q,t;p) maps the 
space of ' BC n - symmetric theta functions of degree m into itself. 

Proof. Every term in the sum is a theta function in each variable with the same functional equation, and thus 
the sum is itself a theta function in each variable. It remains only to show that it is holomorphic, which follows 
from the usual symmetry argument. □ 

Remark. For n = 1, the above operators span a four-dimensional space of difference operators acting on BC\- 
symmetric theta functions of degree m. This 

its action on theta functions |2.'il |2"I] . Moreover, the relations of the Sklyanin algebra are simply given by 
the quasi-commutation relation of Corollary 13.41 below. 

The case m = is worth noting: 
Corollary 1.2. If t n ~ 1 abcd = p, then 

E II d ^Mi0ilMM n ggg n 6(abt^,act^,adt^; P ) (1.20) 

<t£{±1}" l<i<« 6 ^ X i >P' l<i<j<n X j ,P > l<i<n 

Proof. By the proposition, the sum is independent of a^; setting Xi — at n ~ l gives the desired result. □ 
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Remark. This identity is, in fact, equivalent to the identity used in |181 119j to prove BC n -type elliptic hyper- 
geometric identities, a special case of Theorem 5.1 of |29 . . In particular, there is a determinantal proof of this 
identity, which turns out to generalize to a determinantal formula for the difference operator, although we will 
not consider that formula here. 

2 Balanced interpolation polynomials 

For this section, we fix R to be an arbitrary commutative ring. A point on the projective line ¥ 1 (R) is 
represented by relatively prime homogeneous coordinates (z,w) (i.e., the ideal generated by z and w contains 
1), modulo multiplication by units; a polynomial of degree to on P 1 (R) (a section of 0(m)) then corresponds to 
a homogeneous polynomial of degree to in these coordinates. Similarly, given an n-tuple of nonnegative integers, 
one has a notion of a polynomial on P 1 (i?) n of degree (mi,m 2 , ■ ■ ■ ,m n ). The value of such a polynomial at 
a given n-tuple of points on F 1 (R) is only defined modulo units, unless specific homogeneous coordinates are 
chosen for each point; however, this freedom leaves invariant the zero-locus of the polynomial as well as the 
ratio of two such polynomials of the same degree. We will thus feel free to write 

p(xi,x 2 ,...,x n ) (2.1) 

for the value of p at the points Xi,...,x n 6 ¥ 1 (R), with this understanding. 
Of particular importance is the polynomial x ■ y of degree (1,1) defined by 

{z\,Wx)-{Z2,W2)=ZiW2-W\Z2. (2.2) 

Note that x ■ y is invariant under the action of SL/2(i?); since (1, 0) ■ (z, w) — w, we find that x ■ y = if and only 
if x and y represent the same point. Similarly, x ■ y G R* (the unit group of R) if and only if they represent 
distinct points in P 1 (R ) for all homomorphic images i?o of R. 

Definition 3. A symmetric polynomial on V 1 (R) of degree m in n variables is a polynomial on ¥ 1 (R) n of 
degree (m, m, . . . , to) invariant under permutations of the n arguments. 

The space of symmetric polynomials of degree to in n variables (denoted A™) is a free module of rank ( m ^ n ) ; 
this also counts the number of partitions A C to™. An explicit basis of this module is, for instance, given by the 
homogeneous monomials of degree m in the n + 1 multilinear polynomials 

ei(zi, wi, . . . , z n , w„) = ^2 n^n^j) ( 2-3 ) 

/C{l,2,...,n}jSJ j$i 

\i\=i 

< i < n. A polynomial in A™ will be said to be primitive if the ideal generated by its coefficients contains 1; 
note that it suffices to consider its coefficients in the a basis. 
The following decomposition will prove useful in the sequel. 

Lemma 2.1. Fix integers m, n> 0, and a point xq 6 P 1 (i?) (with chosen homogeneous coordinates) . Then we 
have the following short exact sequence: 

(2.4) 



s 



in which the map f is defined by 

f(p)(xi,... ) x n )=p(xi ) ...,x n ) 'Yl {xo-Xi) (2.5) 



Ki<n 



and the map g is defined by 

g(p)(xi, . . . ,x n -i) =p(xo,xi, . . . ,x n -x). (2.6) 
Moreover, this exact sequence splits: there exists a map h : A™_ x — ► A™ such that 50/1 = 1. 

Proof. The claim is clearly invariant under SL2(i?), so we may assume that xq = (0, 1). But then f(p) = e n p, 
giei) = ej for i < n, and g(e n ) — 0; that the sequence is exact follows immediately. The map h is then given by 
h(ei) = a. □ 

Definition 4. Let to, n > be positive integers. A (i?-valued) extended bigrid of shape m™ is a function 

7 : {0, 1} X {1, . . . ,n} X {0, 1, . . . ,m] - P\R). (2.7) 

A bigrid of shape m n is the restriction of an extended bigrid of shape m n to the subset of (a,i,j) such that if 
a — then j < m, and if a = 1, then j > 0. If the extended bigrid 7+ restricts to 7, then j + is said to be an 
extension of 7; in considering the extensions of a given bigrid, we allow the extended values to lie in any ring 
containing R. 

Definition 5. Let 7 be a bigrid of shape to™, and let A C m n be a partition. A balanced interpolation polynomial 
of index A for 7 is a symmetric polynomial P\(', 7) of degree m on P 1 (i?) n such that for all extensions 7 + of 7, 
the following vanishing conditions hold: (1) For all partitions /1 C m™, A (jt /U, 

P* x (... 1 +(0,t^ l )... ]1 )=0, (2.8) 

and (2) for all partitions [i C m™, fj, (f. A, 

P A *(... 7 + (1,*,^:)...;7)=0. (2.9) 

The bigrid 7 is perfect if for all A C to™ , there exists a primitive balanced interpolation polynomial of index A 
for 7. 

In the sequel, we will omit the word "balanced" except when necessary to avoid confusion with Okounkov's 
interpolation polynomials, which we will refer to as "ordinary" interpolation polynomials throughout. (The 
* to denote interpolation polynomials is inherited from shifted Schur functions via Okounkov's interpolation 
polynomials.) The word "balanced" is used for two reasons. First, ordinary interpolation polynomials are 
defined via vanishing conditions together with a triangularity condition; we have replaced the latter by an- 
other, symmetrical, set of vanishing conditions (thus obtaining a more "balanced" definition). Second, the 
corresponding hypergeometric sums (most notably Theorem 14. II below) are multivariate analogues of balanced, 
very-well-poised, hypergeometric series; the corresponding identity for ordinary interpolation polynomials is 
merely very-well-poised. Moreover, although it is easy to degenerate the sums to remove very-well-poisedness, 



9 



there appears to be no straightforward way to obtain a non-balanced sum without degenerating to ordinary 
interpolation polynomials. 

Note that it suffices to consider the generic extension of 7, in which the extended points are all algebraically 
independent over R. Thus the existence of an interpolation polynomial of specified shape is equivalent to 
the existence of a primitive solution to a certain finite system of linear equations, or in other words to the 
simultaneous vanishing of an appropriate collection of determinants. Thus the space of perfect bigrids is a 
closed, i?-rational subscheme of the space of bigrids; moreover, perfection is preserved under base change. 

Proposition 2.2. Let 7 be a perfect bigrid. Then the complementary bigrid 7 defined by 7(0;, = 7(1 — 
a, n + 1 — i, m — j) is perfect. 

Proof. Indeed, an interpolation polynomial of index A for 7 is also an interpolation polynomial of index to™ — A 
for 7, and vice versa. □ 

The simplest example of a perfect bigrid is the following. 

Proposition 2.3. Let 7 be an arbitrary bigrid of shape to 1 . Then 7 is perfect, with associated family of 
(univariate) interpolation polynomials given by 

Pj^l) = II (* -7(0,1,*)) I] (a;- 7(1, l,i)). (2.10) 

0<i<j j<i<m 

Proof. Indeed, the vanishing conditions say precisely that 7(0, is a zero for < i < j and that 7(1, is 
a zero for j < i < m. □ 

This generalizes in several different ways to the multivariate case. The simplest of these is the case of 
"monomial" bigrids, in which 7(0;, i,j) is independent of i. In other words, a monomial bigrid 7 is constructed 
from a bigrid 70 of shape to 1 by j(a, i,j) = 7o(a, 1, j). 

Proposition 2.4. All monomial bigrids are perfect, with associated interpolation polynomials given by 

*--«7)«E II U^o). (2.11) 

7r£S„ l<i<n 

Proof. The given polynomial is clearly symmetric and homogeneous, so it remains to verify the vanishing 
conditions. For this, we observe that if /1 C m n is such that /1 7$ A, then for any permutation tt, [ii < 
for some 1 < i < n. Then the ith factor of the term for it vanishes as required. The vanishing conditions for 
7(1, i, j) follow by the symmetrical argument. □ 

Similarly, the Schur case extends to balanced polynomials. Fix m, n, and let 70 be a bigrid of shape 
(m + n — l) 1 . Then we define a Schur bigrid 7 by: 

l(a,i,j) =7o(a, l,j + n-i). (2.12) 
Proposition 2.5. Schur bigrids are perfect, with associated interpolation polynomials given by the formula 

p *, v _ dct ( jJ A I +»- t (^;7o))i< I j<» det(P;. +n _ i (a;j;7o))i<t,j<n . . 

det(P*_ i (a; j ;7o))i< j , j <„ Ui<i<j< n ( x i ' x i) 
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Proof. We first observe that the numerator vanishes when 

j — ) SO the result is a polynomial; similarly, 
switching Xi and Xj negates both numerator and denominator, so it is a symmetric polynomial. Finally, the 
numerator is homogeneous of degree (m + n— l,m + n— 1, . . . ,m+n— 1), while the denominator is homogeneous 
of degree (n — 1, n — 1, . . . ,n — 1). Thus P£(; 7) as defined is a symmetric polynomial of degree m on P 1 (i?) ra . 

For the vanishing conditions, we observe that, at least in the generic case, the points at which we must 
vanish satisfy x, L ^ xj for all i ^ j. Thus it suffices to show that the numerator vanishes as required; this follows 
by the same argument as in the monomial case. 

Thus Schur bigrids are generically perfect with given interpolation polynomials; the desired result is closed, 
so holds for all Schur bigrids. □ 

Corollary 2.6. Any bigrid of shape l n is perfect. 

Proof. Indeed, any bigrid of shape 1™ is a Schur bigrid. □ 

There is a third important multivariate extension of the univariate case which is not simply an obvious 
generalization of a special case for ordinary interpolation polynomials. In the monomial and Schur cases, the 
polynomials have reasonable formulas, but unlike in the univariate case, do not in general have nice factoriza- 
tions. It turns out that there is a large class of bigrids for which the corresponding interpolation polynomials 
do factor completely. 

Let 77 be a function 77 : {0, 1, . . . , n} x {0, 1, . . . , m — 1} — > P 1 (R), and construct a bigrid 7 by 

7(0,i,j) = ri(i,j) (2-14) 
7(M,j) = r?(*-l,j-l) (2.15) 
A bigrid constructed in this manner will be called a "Cauchy" bigrid. 

Proposition 2.7. Cauchy bigrids are perfect; indeed, we can construct corresponding interpolation polynomials 
via the product expression: 

P A * (*!,...,*„; 7)= 1] II Xi-rtiX'jJ-l). (2.16) 

l<i<n 

Proof. Again, the only nonobvious fact is that P^(;7) as defined satisfies the vanishing conditions. 
Choose fj> C m™ such that fi 7$ A, and consider 

P A *(... 7 +(CM,M,)...;7)= II II 7 + (0,i,^)-r?(A;.,j-l) (2.17) 

l<i<n l<j<m 

Since \l 7$ A, there exists a position 1 < k < n such that fik < Afe. If k' is the largest such k, then we note that 

Afe'+i = /Xfc'+i < fJ-k < A fc < n, (2-18) 

and thus in particular AJ = k for fik < I < Afe. But then 

7+(0, k, n k ) = 7(0, k, Hk) = n{k, Mfe ) = r?(A^ +1 , (/* + 1) - 1), (2.19) 

so the i = k, j = /jfc + 1 factor of the above product vanishes. 

The other set of vanishing conditions follow by symmetry. □ 
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Remark 1. Compare the proof of Lemma 6.3 of |12j. Similarly, interpolation polynomials for Cauchy bigrids 
will be related to the Cauchy identity for interpolation theta functions, Theorem l4. 181 below. In addition, when 
the Cauchy bigrid is also elliptic (see Section EJ, the corresponding interpolation polynomials satisfy a number 
of symmetries, which we will use to prove corresponding symmetries for the general elliptic case. 

Remark 2. This is the only case (at least for a regular (Definition[7J bigrid) for which the associated interpolation 
polynomials factor completely. Indeed, this holds by direct computation for shape l 2 , and every equality required 
of a Cauchy bigrid is supported on some truncation (Definition of shape l 2 . 

The fourth "elementary" case is that of a "delta" bigrid, for which 7(0, — 7(1, whenever both sides 
are defined. If such a bigrid is perfect, then the corresponding interpolation polynomials must satisfy 

PJ(...7 b (i,rt)...;7)K* AM (2.20) 

(here jo(i,j) — 7(0, or 7(1, i,j), whichever is well-defined). Indeed, aside from the fact that the extended 
bigrid 7+ need not satisfy 7 + (0, i,j) — j + (l,i,j), the first set of vanishing conditions says that this vanishes 
unless A C /i, while the second requires that it vanish unless /i C A. That all delta bigrids are perfect requires 
some additional machinery, and is shown in Proposition 12 . 15l below. 

Definition 6. Let 7 be a perfect bigrid. We define four truncated bigrids as follows. The bigrids 7", of shape 
m n , and 7_, of shape (m — 1)", are simply the restrictions of 7 to the appropriate domains; the bigrids ~7 
and _7 are then defined by complementing the restrictions of 7; thus 

-j(a,i,j)=j(a,i + l,j) (2.21) 
_7(a J t,i)=7(a,t > j + l) (2.22) 

The relation of truncation to interpolation polynomials is described by the following propositions. 

Proposition 2.8. Let 7 be a bigrid of shape m", let A C m" _1 be a partition, and let P£ (; 7) be a corresponding 
interpolation polynomial. Then the specialization 

P* x {x 1: x 2 ,..., 7(0, n, 0); 7) (2.23) 

is an interpolation polynomial for 7" of index A. Similarly, if P* ^(j 7) * s an interpolation polynomial for 7 of 
index m ■ X, then 

Pm-xili 1 ' l,m),x 1 ,x 2 , ■ ■ .,x n -x\i) (2.24) 
is an interpolation polynomial for ^7 of index A. 

We cannot quite conclude, however, that if 7 is perfect, then so are 7~ and "7; the constructed interpolation 
polynomials could be trivial even if the original polynomials were primitive. 

Proposition 2.9. Let 7 be a bigrid of shape m n , let A C (m — 1)™ be a partition, and let P^(;7_) be a 
corresponding interpolation polynomial for 7_ . Then 

Y[ x i -7(l,l,m)P;(a ;i ,...,ar n ;7_) (2.25) 

Ki<n 
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is an interpolation polynomial of index A for 7. Similarly, if P^(; -7) is an interpolation polynomial, then 

II aJi- 7(0, »,0)Pi(a!i,...,x n ;_7) (2.26) 

l<i<n 

is an interpolation polynomial of index 1" + A for 7. 

With a slight additional hypothesis, we obtain a converse: 

Proposition 2.10. Let 7 be a bigrid of shape to™, and let A C (to — 1)™. If the point 7(1, 1, to) has unit inner 
product with all points in the image of 7_ , then all interpolation polynomials of index A for 7 are obtained 
from interpolation polynomials for 7_ via the above proposition. Similarly, if the point 7(0, n, 0) has unit inner 
product with all points in the image of _7, then all interpolation polynomials of index 1™ + A for 7 are obtained 
via the above proposition. 

Proof. Passing to the complementary bigrid reduces the second claim to the first. Thus consider an interpolation 
polynomial P£(;7). The partition m is not contained in A, and thus the polynomial 

P A *(7(l,l,m),x 2 ,...,x„;7) = 0. (2.27) 

It follows that P£ (xi , . . . , x n ; 7) is a multiple of x\ ■ 7(1, 1, to), and thus, by symmetry, that it is a multiple of 

II asi^M.m). (2.28) 

l<i<n 

We need to show that 

Pxh(^A,m),x 2 ,...,x n ;-f) 

Ul<i<n X i -7(1,1,™) 

satisfies the vanishing identities required of an interpolation polynomial for 7_; by the hypothesis, it suffices to 
show that P^ (; 7) itself satisfies the identities. But the required identities are a subset of the vanishing identities 
that P^ already satisfies. □ 

Corollary 2.11. Let 7 be a perfect grid of shape to™. 1/7(1, 1,to) has unit inner product with the image of 
7_, then 7_ is perfect; 1/7(0,71,0) has unit inner product with the image of _7, then _7 is perfect. 

Naturally, the interesting cases of perfect bigrids are those for which the interpolation polynomials are 
uniquely determined. This is an open condition (the complement of a closed subscheme again cut out by 
determinants); unfortunately, the ideal for the complement appears to be prohibitively complicated. We thus 
need to find a simpler set of sufficient conditions for uniqueness. As in ^3 Definition 2.1], the idea is to choose a 
subset of the equations and add a normalization condition in order to obtain a square system of linear equations 
with a simple determinant. 

Definition 7. A bigrid is regular if for all i > i' , j < f, a, [3 S {0, 1}, 

"f(a,i,j)-j(/3,i',j')&R*. (2.30) 
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Note that regularity of a bigrid 7 implies that the extra hypotheses of Proposition ^. lOl hold for all truncations 
of 7; regularity is slightly stronger in that the additional conditions 

7 (M,j)- 7 ((M',j-l)G J R* (2.31) 

for i > i' must also hold. (It is unclear whether these extra conditions are necessary for primitive interpolation 
polynomials to be unique and form a basis; they are needed for our proof below, however.) It follows therefore 
that if 7 is a regular perfect bigrid, then so are 7_ and _<y. 

Given a bigrid 7 of shape to™ (to > 0), a partition A C to™, and an integer < I < n, we dehne a point 
71(A) eP 1 ^)" by 

7l(A) = ( 7 ([1 < q,l,Ai),7([2 < Z],2,A 2 ),...,7([b < /],n,A„)), (2.32) 

where [i < I] is 1 if i < I and otherwise. For a pair of partitions A,/i C to™, we define jx(p-) = 7z(/-t) with Z 
determined as follows. Let Zo be the largest index such that m ^ A; (Zo = if no such index exists), and let l\ 
similarly be the largest index such that fii = to. If Iq = or fj,i < A/ then take I = l\\ otherwise, take I = Iq. 
Thus, for instance, for to = n — 3, A = 21, we have: 



721(0) = 


(7(0, 1,0), 7 (o, 2,0; 


,7(0,3,0)) 


(2.33) 


721(1) = 


(7(0,1,1), 7(0,2,0; 


,7(0,3,0)) 


(2.34) 


721 (2) = 


(7(1,1, 2), 7(0,2,0; 


,7(0,3,0)) 


(2.35) 


721(H) = 


(7(0, 1,1), 7 (0, 2,i; 


,7(0,3,0)) 


(2.36) 


72i(21) = 


(7(1, 1,2), 7(0,2,i; 


,7(0,3,0)) 


(2.37) 


72i(22) - 


(7(1, 1,2), 7 (1, 2, 2; 


,7(0,3,0)) 


(2.38) 


72i(Hl) = 


(7(1,1,1), 7(1, 2,r 


,7(1,3,1)) 


(2.39) 


72i(211) = 


(7(1, 1,2), 7(1, 2,i; 


,7(1,3,1)) 


(2.40) 


72i(221) = 


( 7 (1, 1,2), 7(1, 2, 2; 


,7(1,3,1)) 


(2.41) 


72i(222) = 


(7(1, 1,2), 7(1, 2, 2; 


,7(1,3,2)) 


(2.42) 



Definition 8. Fix a ring R, positive integers to, n, a R- valued bigrid 7 of shape to™, and a partition A C to™. 
A primitive symmetric polynomial p of degree to and n variables is a quasi-interpolation polynomial of index A 
and bigrid 7 if for all partitions fi C to™, fi =/= A, 

p(7A(A0) = °- ( 2 - 43 ) 

Lemma 2.12. If the bigrid 7 of shape to™ is regular, then for all A C to", there exists a quasi-interpolation 
polynomial of index A, unique up to scale. The quasi-interpolation polynomials for 7 span the space of symmetric 
polynomials of degree to on P 1 (i?)™. 

Proof. Introduce a total ordering <; on partitions by first sorting the parts into increasing order, then lexico- 
graphically ordering the sorted partitions. We construct a corresponding filtration of A™ as follows. (That is, a 
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sequence of submodules F fJi (X;'y) such that F M (A;7) D F^(A;7) whenever /i <; f.) Recall from Lemma r2.il the 
short exact sequence: 

O-A^-^A^ASU-O, (2.44) 

in which the map / is defined by 

/(p)(xi,...,x„) = p(xi,...,x„) JJ (7(0, n, 0) • x 4 ) (2.45) 

l<i<n 

and the map g is defined by 

g(p)(xi, x n _i) = p(x x , . . . , x n _i, 7(0, n, 0)). (2.46) 

The submodules F^X;^) c A™ are then defined as follows. First, as a base case, i*o(A;7) is always taken to 
be A™; this in particular covers the case m = or n — 0. If /i„ = A„ = 0, rn > 0, then 

F M (A; 7 )= ff - 1 (F At (A; 7 -)); (2.47) 

if /Lt„ = but A n > 0, then 

F^Xll) = g-^F^')), (2.48) 
where 7'(a, i, j) := 7(0, z, j) is a delta bigrid of shape m" -1 . If £i n > 0, A„ > 0, then 

F /1 (A; 7 ) = /(i ?, /1 -i«(A-l";-7))- (2-49) 

Finally, if /i„ > but A n = 0, then 

^(A; 7 ) = /(^-i"(0;7")), (2-50) 

where 7" (a, i,j) := 7(1,2, j — 1) is a delta bigrid of shape (m — 1)". 

Moreover, since the short exact sequences used to construct the filtration are all split, we can construct a 
corresponding sequence of polynomials p M (; A; 7) such that 

F fl = (p u : n<iv). (2.51) 

(When m = or n = 0, we take po = 1) From the construction of this filtration, it follows easily that 

P m ( 7 aM;A; 7 ) =0 (2.52) 

if f <; fi, while 

P M ( 7 a M ; A; 7) eiT, (2.53) 

using the fact that 7 is regular. (Note that if 7 is regular, then so are the delta bigrids 7', 7" used above) In 
particular, the determinant of the matrix 

(p/.(7A(f))Wcm» (2-54) 
is a unit, so for any collection of values € R, there is a unique polynomial p such that 

?(7a(m)) = / M . (2-55) 
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Taking f^ = Sx^ gives the desired existence and uniqueness claims for the quasi-interpolation polynomials. 
For v <i A, "fx(v) is independent of A; it follows that 

^(7m"M)=0 (2.56) 

when v <; A, while 

Px(lmn(X)) €R*. (2.57) 
Thus the corresponding matrix has unit determinant, which immediately implies that the P£ form a basis of 

K- □ 

Remark. Compare the proof of Proposition 2.6 of |13|. 

Theorem 2.13. Let 7 be a regular, perfect bigrid of shape m n . Then the corresponding primitive interpolation 
polynomials are uniquely determined up to scale, and form a basis of the space of symmetric polynomials of 
degree m on P 1 (i?)". Moreover, for all \, 

P A *(... 7 + (0,i,A0...;7) (2-58) 
is a primitive polynomial in the indeterminates 7 + (0,z,m) ; 1 < i < n. 

Proof. We find that the vanishing conditions defining the quasi-interpolation polynomials form a subset of the 
conditions defining the interpolation polynomials, so each interpolation polynomial must be a unit multiple of 
the corresponding quasi-interpolation polynomial. We then find that 

P x *(... 1 + (0,i,\ i )...; 1 )eR* (2.59) 

in the special case 

j + (0,i,m) = 7(1, i,m), (2.60) 

which implies primitivity in general. □ 

Corollary 2.14. If 7 is a regular, perfect bigrid, then the same is true for any truncation 0/7. 

Proof. The only thing remaining to prove is that if A„ = and P£ (; 7) is a primitive interpolation polynomial 
for 7, then substituting 7(0, n, 0) for one argument gives another primitive polynomial. But the corresponding 
statement is true for quasi-interpolation polynomials. □ 

In the sequel, we will restrict our attention to regular, perfect bigrids. 

Proposition 2.15. Delta bigrids are perfect. 

Proof. Since regular delta bigrids are Zariski-dense in the space of all delta bigrids, it suffices to consider the 
case of a regular delta bigrid 7. Let -P\(;7) denote the corresponding quasi-interpolation polynomials; thus 

P A *(...7(0,i,/i*).-.;7) = <V (2.61) 
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Taking a ring extension as necessary, we can also choose a regular Cauchy bigrid 7' such that j'(0,i,j) = 
7(0, with corresponding interpolation polynomials. But then 

P;(... 7 + (0,^)---;7') = (2.62) 

unless XC.fi; thus the polynomials corresponding to 7 are triangular with respect to the basis P£(;j') under 
the inclusion ordering. In particular, 

^a(;7') = E c V^(;7) (2.63) 
for appropriate coefficients ca^, and therefore 

^(;7) = E d v^(;7') (2-64) 

^DA 

for another set of coefficients dx^. But then if fi 2 A, 

P A *(... 7 + (0, i ,^)---;7)= E ^P;(...7 + (0,i,/i i )---;7) = 0. (2.65) 

The other set of vanishing conditions follow symmetrically; the result follows. □ 

Given a pair (C, cj>), where C is a curve of genus 1 and is a quadratic map C — > P 1 , the line bundle 0(1) 
on P 1 pulls back to a line bundle on C, on which the Galois group Z/2Z acts. Thus symmetric polynomials 
on P 1 pull back to sections of the corresponding line bundle on C™, invariant under the action of BC n . When 
C is a complex elliptic curve and preserves the identity, symmetric polynomials of degree m in n variables 
thus pull back to i3C n -symmetric theta functions of degree to. In the next section, we will construct certain 
"interpolation theta functions", pull-backs of interpolation polynomials associated to "elliptic" bigrids. We will 
consider the case of general curves C in Section below. 



3 Interpolation theta functions 

Definition 9. Fix integers m, n > 0, a partition A C m", and generic points a, b,q,t € C*; also choose another 
generic point »eC. The interpolation theta function 

P* {m ' n \ Xl , . . .,x n ;a,b;q,t;p) (3.1) 

is the unique function on (C*)" x (C*) 4 with the following properties: 

(1) P\ ™'" is a BC„-symmetric theta function of degree to in the variables X\ , . . . , Xn • 

(2) For any partition fi c to™, fi 7^ A, let lg be the largest index such that fii 7^ Az, let Zi be the largest index 
such that fii — m (0 if none exists). If fii a < A; , then set I = l\, otherwise set I = Iq. Then 

p*(m,n) (&(fl - w ; ^ bqm -„ lt l-l jaqfil+H n-l-^ _ ^ fc; ^ f;p) = q (3.3) 
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(3) p*( m '"' is normalized by the specialization 

P* (m ' n) (. . . vt n -> ...;a,b;q,t;p) = C^f^av, a/v; q, trfC&^^bu, b/v; q, t;p) (3.3) 

Indeed, aside from the normalization, this is simply the pull-back of the definition of quasi-interpolation 
polynomials relative to the bigrid 

7 (0, := HaqH^) ^(bq^f- 1 ), (3.4) 

where <fi maps the elliptic curve C* / (p) to P : (C) via a pair of BCi-symmetric theta functions of degree 1; i.e., 
for some c, d € C* , 

d>(z) = (e(cz ±1 ;p),e(dz ±1 ]P )). (3.5) 

Thus by Lemma ?2. 121 interpolation theta functions exist and are unique, at least for generic values of a, b, q, t, v. 
We note in particular that the associated bigrid is regular precisely when 

q j t i iabq m - 1 t n - 1 (p) 0<i<n,0<j <m (3.6) 

(b/a)qH l i abq m -H n - l {p) \i\<n- 1, |j| < m - 2 (3.7) 

q J t l <£ (p) < i < n, 1 < j < m (3.8) 

(b/a)qH l i (p) \i\<n- 1, |j| < m - 1 (3.9) 

q j t i i {ab)- 1 q 1 - m t 1 - n {p) <i <n,2< j <m (3.10) 

(b/a)qH l £ (a6)- 1 g 1 - m t 1 - n (p> |t| < n - 1, |j| < m - 2 (3.11) 

Moreover, we can also conclude that p^ m '")(; flj 5; ^ i s a meromorphic theta function in each of a, 6, g, t, v. 
(We will show below that it is independent of v, explaining its omission from the notation) 

As one might expect from the fact that we are considering these functions at all, interpolation theta functions 
indeed satisfy extra vanishing conditions corresponding to the fact that the associated bigrid is perfect. The 
normalization can thus be justified as follows: viewed as a function of v, the left-hand side must (by extra 
vanishing) vanish wherever the right-hand side vanishes. By degree considerations and the theta function 
condition, it follows that the ratio must be independent of v. 

As in we first show that the interpolation theta functions satisfy a difference equation, and in the 
process prove independence of v. 

Theorem 3.1. Let c, d £ C* be chosen so that q m t n ~ 1 abcd — p. Then 

(B) (a, b, c, d; q, t; p)P* x {m ' n) (; ^qa, ^db; q, t; p) 

= Y[ e(abq m t n -\acq x 't n - l ,bcq m - x 't l - 1 ;p)Pl {m - n \;a,b;q,t;p) (3.12) 

l<i<n 

Proof. To begin with, we take the v normalizing the interpolation theta function in the left-hand side of the 
equation to be ^fq times the v normalizing the interpolation theta function on the right. 
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We first need to show that the left-hand side satisfies the relevant vanishing conditions. Suppose more 
generally that we evaluate the left-hand side at a point of the form 

bq m -^,. . . , bq m -> J -n l - 1 ,aqM+H n - 1 - 1 , aq^ (3.13) 

for an arbitrary choice of Z, where /x simply satisfies the conditions 

m > Mi > M2 > •• • > Hh (3-14) 
W+i > > • • • > A*n > 0, (3.15) 

but need not satisfy fil > m+\. Carrying this point through the difference operator and using the BC n symmetry 
of P*( m ' n > : we find that in each term, p*( m >")(; ^Jqa, *Jqb; q, t;p) is evaluated at a point of the same form, with 
a,b replaced by <? 1 ^ 2 a, q 1 / 2 b, and p replaced by v satisfying 

tH < Vi < m + 1, 1 < i < I (3.16) 
Pi — 1 < Vi < /ii, I + 1 < i < n (3-17) 

(Indeed, i>i = \ii + 1/2 — <J% if 1 < i < I, and Vi = fa — 1/2 — <Ti if £ + 1 < i < n.) We furthermore observe that 
if m > i/i, < Vi+\ for z 7^ £, or v n < 0, then that term of the expansion necessarily vanishes. Thus the only 
surviving v satisfy 

m > v x > i/ 2 > ■■■ > v u (3.18) 

n+i > n+i > ■ ■ • > v n > 0. (3.19) 

In particular, if (i and I are chosen to make the right-hand side vanish, then all terms on the left-hand side will 
also vanish, as required. 
It follows that 

£>W (a, b, c, d; q, t ; p)P* x M (; y/qa, ^qb; q, t; p) oc P^ m > n \ ; a, b; q, t; p), (3.20) 

for some scale factor independent of Xi, . . . , x n . Taking p = A, and I to be the largest index such that A; = to, 
we find that the interpolation theta function on the right has a nonzero value (by Lemma l2.12[l . and that only 
one term on the left survives. The desired dependence of the scale factor on c follows immediately. 

If we now take c = v and set Xi = vt n ~ l , again only one term survives on the left, and the difference equation 
follows. By symmetry, the same difference equation would hold if we had normalized the interpolation theta 
function on the left using q~ 1 / 2 v instead, and thus these interpolation theta functions agree. In other words, 
interpolation theta functions are invariant under v i— > qv; since q is generic, they must indeed be independent 
oft). □ 

Remark. In particular, we find that the interpolation theta functions are solutions to a "generalized eigenvalue 
problem" (ala [H]), to wit: 

^ (a ' 6 ' C ' q^/^abc 5 q ' <; P) P A (m ' n) (; V^a, y/qb\ q,t;p) (3.21) 



, Ki<n 



jj eiacq^t^Mq™-^- 1 ^) \ (n) p , p *(m,n), r rh , s 

11 , { adq^M' q -^-p) ) D {aAc > q^abc> ;q ' t;P)P > ^^^^Hp) 
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Corollary 3.2. Fix a partition X C m n , and let 1 < / < n and a sequence < [i n < < • • • < m+i be 

chosen such that fii+i < Xi+i- Then 

P* (m ' n) (xi, . . . , x u aq^H n ~ l ~\ aq^;a, b; q, t; p) = 0. (3.22) 

Similarly, if m > /zi > fi2 > ■ ■ ■ > (J-i is chosen such that fj,i > A; , then 

p*(m,n) (6gm - Ml ^ bq m-^ t l-l :Xl+u . . . , ^ a, fe; g? t;p ) = . (3.23) 

Proof. By analytic continuation, it suffices to prove that if /ii, .. . ,fi n are as above, such that either \xi > Xi or 
fJ-i+i < Xi+i, then 

p*(m,n)( bq m-ni ^ ^ bq m-^ f l-l ( aq n l+H n-l-l > ^ ^ . ^ & . ^ t;p j = q (3 24 ) 

We proceed by induction on 

]T(m-^)+ £ w . (3.25) 

l<i<Z 2 + l<i<ra 

Indeed, if we evaluate both sides of the difference equation l|3.12|l at this point, we find by the inductive 
assumption that at most one term on the left survives (that for which 0$ = 1, 1 < i < n). If > A;, 
then setting c = bq" 1 ^^ 1 ^ 1 makes this term vanish without annihilating the factor on the right; similarly, if 
M+i < A; + i, we set c = aq t ^ l+1 t n ^ 1 ^ 1 . The result follows. □ 

Corollary 3.3. We have the symmetry property 

P^(; a, b; q, t;p) = P* X M (; b, a; q, t;p) (3.26) 

Proof. This follows immediately from Corollary 13.21 and the fact that the normalization has this symmetry. □ 

This also gives the following commutation relation for the difference operators: 

Corollary 3.4. // cd = c'd' , then 

D {n \a, b, c\ d'; q, t;p)D™ \q l ' 2 a, q 1/2 b,q" 1/2 c, q^ 2 d; q, t;p) (3.27) 

= r»(a, b, c, d; q, t;p)D™ {q l ' 2 a, q^ 2 b, q~V*d , q^d 1 ; q, t;p) 

Proof. Suppose q' m t n ~ 1 abcd = p. Applying both sides to interpolation theta functions, we find that the relation 
holds on the full space of _BC„-symmetric theta functions of degree m. Now, applied to a _BC n -symmetric 
meromorphic function /, both sides can be expressed in the form 

J2 C e (...x i ...;a,b,c,d,c',d')f(...q ei x i ...) (3.28) 

ee{-l,0,l} n 

for suitable functions C e ; the claim is that the two resulting families of functions agree. For m sufficiently large, 
this follows from the action on theta functions; since each C e is a meromorphic theta function, the result follows 
by analytic continuation. □ 
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When q' m t n ab = pq, the associated bigrid is a Cauchy bigrid; we thus obtain the following formula for 
interpolation theta functions in that case: 

Proposition 3.5. If q rn t n ab = pq, then 

P* x {m ' n \...x t ...;a,b;q,t;p) = J[ 9(at n -^ q^x^at^ q^ 1 / 'x t ;p), (3.29) 

l<z<n,l<j<m 

We will be using this Cauchy specialization quite frequently; the point is that by using the difference 
equation, one can shift ab arbitrarily, but cannot change a/6. Thus the Cauchy specialization, which includes 
cases having all values of a/6, gives a suitable base case for inductive arguments. (It thus plays much the same 
role that Macdonald polynomials did in [T?)].) For instance: 

Proposition 3.6. Interpolation theta functions satisfy the identity 

P* im ' n \...aq x H n - i ...;a,b ] q,t ]P ) (3.30) 



w> , PQ , ^+ ( tn ~ lb * ^{pq,t-q,t-p)C+(t 2n - 2 a 2 -q 1 t-p) 
= C ra--x{-^^,t-p)C+ n _ x {——- ) q,t-,p) T^TI^Z Z3 



^Of,'*' m»-AV q m a >*> "J C° x (t n ;q,t;p) 

Proof. Using the difference equation, it follows immediately that this formula holds, up to a scale factor de- 
pending only on a/6 (more precisely, independent under (a, 6) i— > {^/qa, ^/qb), but by analytic continuation, this 
amounts to the same thing); taking the Cauchy specialization shows that this scale factor is 1. □ 

The other important base case is the delta case: 

Proposition 3.7. Ifq m t n - 1 ab= 1, then 

P* x {m - n \...aq^t"-\..;a,b;q,t;p) 

-S,C° (vat'^-a tv)C+ ( 1 . a t v] C x (P9,t;q,t;p)C+(t 2n - 2 a 2 ;q,t;p) 
-^C mn _ x (j,qt >q ,t,p)C mn _ x (-^ >q ,t,p) C° x (t^q,t;p) 

The dependence of interpolation theta functions on m is fairly simple; we have the following identity. 

Proposition 3.8. If X C m", then 

P x m+k ' n \-.-x t ...;a,b;q,t;p)= JJ_ 0{bxi,b/xi;q;p) k P x ^ m ' n \- ■ -x t . . .;a,bq k ;q,t;p). (3.31) 

l<j<n 

In particular, 



P^ rn ' n '(...x i ...;a,b]q,t]p)= Yl 0(bx i ,b/x i ;q- 1 p) m (3.32) 



* (m.n) , 

. . . a,, . . . , u, u, q, u,pj — 

l<i<n 

In particular, the dependence on m is independent of A, and thus the following abelian functions are well- 
defined: 

Definition 10. The interpolation abelian function R* x n \;a,b;q,t]p) is defined by 

P^ m ' n \;a,q-^b;q,t;p) 

for any m > Ai. 
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Of course, the price we pay for dealing with these functions is that we cannot evaluate them at points 
bq tLi t n ~ t 1 and as a result the complementation symmetry becomes more complicated. However, the abelianness, 
the lack of dependence on to, and the resulting ability to have an infinite family of functions, is often an overall 
win, and as a result we will tend to use these functions in preference to the interpolation theta functions. 

Translating the above identities to the interpolation abelian functions, we obtain: 

Proposition 3.9. If t n ~ 1 abcd = p, then 

D {n) (a, b, c, d; q, t; p)R x {n) (; y/qa, y/qb; q, t; p) 

= Yl e{abt n -\acq X H n -\bcq- X H l - 1 -p)Rl {n) {-a,b-q,t-p) (3.34) 

l<i<n 

As an abelian function in v, 

R< n \ v t n -* ;a ,b;q,t;p) = A^" V&l*" _1 <™, <V«; <?> (3.35) 

Ift n ab — pq, then 

p .(n), . . x TT ^"-V-^ar-V-V^ip) 

R x y >{x;a,b;q,t;p) = — — r (3.36) 

J -- L 6{at n q3 L Xi,at n qJ L /x l ;p) 

Finally, 

We also note the following symmetries that follow from symmetries of the definition of interpolation theta 
functions. 

Proposition 3.10. The interpolation abelian functions satisfy the following symmetry identities: 

Rt\- ..z i ...;a,b;q,t;p) = ( Q fc2 )M t -*»W^ )#*(«)(. l/ a , 1/6; 1/q, l/t;p) 

(3.38) 

R* x {n) (. . . Zi . . . ; a, b; q, t; p) = R* x [n) { Zi . . . ; -a, -b; q, t; p) (3.39) 

R* x {n \. ..y/p Zi ...;y/pa, y/pb; q, t;p) = (bt^/aq^t^q-^R^i. . . z t . . . ; a, b; q, t;p) (3.40) 
R*\ n) (- ■ ■ Vpzi • • • ; \/pa, b/y/p; q, t;p) = {pq/ab) W R* x {n) (. . . z t . . . ; a, b; q, t;p) (3.41) 

R*^ +X (...z i ..,,a,b;q,t;p)= J] gi R* x ^ (. . . z, . . . ; of", 6/<T; q, t; p) (3.42) 

l<i<n a \ b Z i ' 1>P)m 

R-r% ..*...,„,„<,..., „<*-;«, M. «: P) - 7 ^ i ^0- )R TH. t-u, * ,, * p) (3.43, 

4 Binomial coefficients and hypergeometric identities 

Definition 11. The generalized binomial coefficients are meromorphic functions of a, b, q, t defined by 

( X ) ■■= A„(£|t n , l/b;q,t;p)ltf% . . . . .-^^b/y/E^^p), (4.1) 

V/V [a,b];q,t;p 
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for any integer n > £(X), t{p)- 

Note that by equation i|3.39|l this is indeed independent of the choice of i/a. While for many purposes the 
given normalization of the binomial coefficients is the nicest, it does have the significant drawback of being 
singular at points of the form b = q~ k t l for integers k, I > 0. With this in mind, we introduce a second 
normalization: 

'A\ _ Al(a\b,v 1 ,...,v k ;q,t;p) fX\ 



\/V [ a ,b](v 1 ,...,v h ); q ,t; P ' A ° ia/b\l/b, ,v k ;q,t;p) V/'/ 

(We will omit the parentheses when k = 0.) Thus for instance for b = 1 we find 

'A 



= <V (4.3) 

[a,,i\(vi;--.,v k );q,t;p 



The standard normalization satisfies particularly nice transformation laws: 

'A\ /A N 



Av [pa,b];q,t;p vAv [a,b];q,t;p 

") = Q 

Av [a,pt];gr,t;p VAV [o,6];g,i;p 

A\ / A N 



vAV [l/o,l/6];l/a,l/i;p \AV [a,6];g,t;p 

We also note the following special values: 

A 




(4.4) 
(4.5) 
(4.6) 



= 1 (4-7) 

[a.b]:q,t;p 

A\ C+(a;q,t;p)A° x (a/b\l/b;q,t; P ) 



A/ [a,b];q,t; P <?+ (§; g, t; p)A° A (a|6; q,t; p) 

ffl 71 \ n 

= A x (-\f\q-"\t 1 - n q m a,l/b; q ,t;p) (4.9) 
A / [o,6];g,i;p 

Furthermore, we readily obtain the following symmetry, from the corresponding symmetry of interpolation 
functions: 

l™"+<J[a,6 ]; ,, t;P = A° (q 2m a\b, pqaq m /b, pqt^q™ , t 1 - V m a; g, t; p) /A\ 
(m»)[a,6] ;g ,t ; p A°(g^o/6|l/6, W a<r,l»8t^ 1 8 m ,t 1 - B g 3m a/6;? ) t;p) W [g ^a, 6 ] ; ,, t;P 

As announced in |15) . the main summation identity for elliptic generalized binomial coefficients is the 
following identity; as the proof below suggests, it can be thought of as a "bulk" difference equation, ala the 
bulk Pieri identities and branching rule of [J6J. 

Theorem 4.1. For any partitions k C A, and generic parameters a, 6, c, d, e, g, t S C* smc/i i/iai frccfe = apg, 

= A0 / i U h 71 r~T~ 1^ & (a/b\c/b,pqa,d,e;q,t;p)[ 

K J [a,c]; 9 ,t;p A^(a|c, »d, be,pqa/b; g, t;p) ^c^tcA W [a,b] ;g ,t ; p W [ a /6, c /6] ;gA p 

(4.H) 
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Proof. If we take the difference equation (|3.12l) for interpolation abelian functions and rewrite it in terms of 
binomial coefficients, we obtain an identity of the form: 

'A\ A a R {aJb\lJ ab, ac, ad,pq 2 a 2 ; q, t;p) 



(4 12) 

K J[a^ab];q,t; P A°(a 2 |afe, ac, ad; q, t; p) 



^2 A ° (i a2 \<i ab , 9 ac > i ad \ iA p) f x/h (a;q,t;p) 



^ ' [a 2 q,abq];q,t;p 

with abed — p; the coefheients F^/^ia; q, t;p) are complicated, but are in principle explicit. Extending F to be 
if \x -fc A, we can replace the range of summation with k C [i C A without changing the sum. Now, in the special 
case b — 1/aq (and thus cd — pq), the binomial coefficient on the right-hand side becomes a delta function, and 
we thus obtain the identity 

F x/K .(a;q,t;p) = ( X ) . (4.13) 

W [a 2 ,l/q];q,t;p 

Substituting this in gives the special case b = 1/q of the theorem. 

We next observe that the values of b for which the theorem holds is closed under multiplication. Indeed, if 
b = b\b 2 such that the theorem holds for b = b\ and b = b 2 , then wc in particular find using the case b = b\ of 
the theorem that 

'A\ A°(a/bib 2 \l/bib 2 ,pqa/bi,bid, ce; q, t; p) 



- (4-14) 
M/ \aMi*];q,t; P A° x (a\bib 2 ,pqa/bi,bid,ce;q,t;p) 



^2 A°( a / 6 il 6 2,P3M,ce/6i;g,i;p) 



A ) C) 

fV [o,6i];g,t;p V K / [a/bi,6 2 ]; 9 ,*;p 



Substituting this in to the right-hand side, we find that the resulting double sum can be summed over /1 using 
the case b = b 2 of the theorem, and then summed over v using the case b = b\ again. 

We thus conclude that the theorem holds whenever b = q~ k for k a positive integer. Since for a, c, d fixed, 
both sides are meromorphic theta functions in b with the same multiplier, the result for general b follows by 
analytic continuation. □ 

Remark 1. This identity includes a number of known elliptic hypergeometric identities as special cases. Taking 
both A and k to consist of a single part gives the Frenkel-Turaev summation for univariate elliptic hypergeometric 
series [§] (the elliptic analogue of the s4>7 summation due to Jackson). Taking A to be a rectangle and n = 
gives an identity conjectured by Warnaar fZ§\ Corollary 6.2] (also proved in US]), which in our notation reads 

a 1 un -m l l u h + \ Cm- (W«, pqa/bpbi , pqa/b b 2 ,pqa/bib 2 ; q, t; p) 

V A M (at ,q ,b ,b 1 ,b 2 ,b 3 ;q,t;p) = —f- — — — , (4.15) 

L-m" {pqa/bo,pqa/bi , pqa/b 2 ,pqa/b b 1 b 2 ; q, t; p) 

where q~ m t n bobib 2 b 3 = a 2 pqt. 

Remark 2. When t — q, the binomial coefficient is related to Schur-type interpolation functions, and can thus 
be expressed as a determinant; to be precise, we have: 

A\ = n dct ) {m 

\M/r m 11 q- x <9(q x *- i - X 3+3,q x i+ X J+ 2 - i ~3a;p) Ki,j<n \ + n-J/< a ,1 / 

\f"/ [a,b];q,q;p l<i<j<« H Vy ' y ' Fl ~ \ ■> I [ 2 °_ 2 ,b];q,q;p I 
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(The given scale factors can be derived via Warnaar's determinant identity, Lemma 5.3 of |55]; also note that 
setting X — fi makes the determinant triangular, so trivial to evaluate.) This gives rise to Warnaar's Schlosser- 
type identity (Theorem 5.1 of |29|). via the case t = q, b = q~ k , K = above, in the following way. When 
thus specialized, the binomial coefficients with lower index k are ratios of C° symbols, while the remaining 
binomial coefficient can be expressed as a determinant. Now, for this coefficient not to vanish, we must have 
Aj — k < Hi < Aj for each i; if in addition A; — A,; + i > fc, the relevant matrix becomes diagonal. Since the sum 
is terminating for generic A, we can analytically continue in q Xi , 1 < i < n (summing over A, — /ij); the desired 
identity results. 

Corollary 4.2. The binomial coefficient ( x ) vanishes unless u -< A, in which case 

W la,l/q];q,t;p ^ 

A\ yr 8(q^-H 2 -^- l a-p) yr 9 (q x ^ - 1 pq; p) 

N { a,l/ q] ; q ,t; P ^ e(q^-H^-H;p) ^Oiq^-H'-^apqy^p) 

(i J) 6 , ! : ' 2 " '">r-P! 9{q^-H^pq;p) 

Proof. The proof of Theorem 14.11 gives a formula for the binomial coefficient, which in particular implies that 
it vanishes when required. Simplifying this as in the proof of Theorem 4.12 of ^3] gives the desired result. □ 

Remark. Similarly, for any integer k > 0. ( A ) r vanishes unless a -<k A. 

N A t/ [a,q~ k ];q.t\p 

When c = 1 in the bulk difference equation (|4.11|l , the left-hand side becomes a delta function, and we thus 
conclude: 

Corollary 4.3. Elliptic generalized binomial coefficients satisfy the inversion identity 

KCMCA W [a,b];q-t;p W [a/b,l/b];q,t;p 

From the form of Cauchy interpolation polynomials, the following symmetry of interpolation theta functions 
is essentially automatic: 

P;(™'"> (A; a,pq/q m t n a- q, t; p) = (n m - //; ^p/a, ^pf 1 ' 1 q m a; l/t, 1/q- p) ; (4.19) 

note in particular that the roles of argument and index are interchanged via this symmetry. The special form 
of the bigrid in our case gives rise to another symmetry as well, which follows from a computation along the 
lines of Lemma 2.1 of |16j : 

P*(m,n) (x . n nn/n m f n n . n f . n) _ C° „ ft"; g, t;p) C^f^O*, (pg/t)f"; q, t;p) Cgft""^ 2 , g~"; g, t;p) 
„ [A^a^pq/q ,q, ,p, ^^.^ C°(t^ , (pq/t)q-™; q,t; p) C^t^a 2 ,t-;q,t; P ) 

^(».™)( A ' ;tn -i/2 ? m-i/2 ajPg 8/a t i/a /a . 1/fj l/q . p) (4 2Q) 

We thus obtain an action of the group Z\ on Cauchy interpolation theta functions; this symmetry, it turns out, 
actually extends to all interpolation theta functions. In terms of binomial coefficients, we obtain the following 
two symmetries: "duality" and "complementation" . 
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Corollary 4.4. (Duality) 

'A\ / A' 



, , , , (4- 21 ) 

fV [a,b];q,t;p / [aqt,b];l/t,l/q:p 

Proof. When b — pq/t, the relevant interpolation functions are of Cauchy type. We thus find that 
'A 



= F\(a;q,t;p)G t _ t (a;q,t;p)( J 

,t;p V* / fan 



/V [o,P9/t];«,t;p V'- 1 / [ogt,pg/t];l/t,l/g;p 

for appropriate functions F\, G M , with Go = 1. Setting jj, = shows that Fa = 1; setting /x = A then shows 
that G M = 1. (These could, of course, also be computed directly.) 

The bulk difference equation l|4.11|) is preserved by this symmetry; it follows that the set of b for which this 
duality holds is closed under multiplication, and thus the result follows by analytic continuation. □ 

Corollary 4.5. The binomial coefficient ( x ) , , vanishes unless a -<' A, in which case 

^ [a.£];g,£;p 

A\ n 9(g^-H 2 - x '^a:p) yr g^g^Hjjg) 

M[a,t]„*P { £)\x 0{q^-n^- l pq-p) ^ 9^+^-^ apq/t*; p) 

yr 8(q x >-n x 'i-*pq;p) yr 9 (q x < +^ t 2 - X 'i -' pqa / f; p) 
e(q^-H^a/t;p) ^ e(q^-H^H;p) 

Remark. Similarly, ( x ), , vanishes unless a <' h A. 

J ' W [a,t k ];q,t;p ^ h 

Corollary 4.6. If Ai,/Zi < m, the binomial coefficients satisfy the identity 

( m l' X ) A°(-g-\h 22£ Pi 1 m °- ■ n f-n) /\\ 

(m»)[a,6];9,t;p g™t"+i ' t^™-"b ! 9, P) VV [t-2»a,6]; 9 ,t;p 

Similarly, 
Corollary 4.7. (Complementation) 

G)[q,6] ;g ,t ;P = A M (q/fc|t", <T"', 1/6, g m t 1 -"a; g, t;p) An" - /A 

A m . (a|t n , <r m , 6, g m t 1 -"a/6; ?, t;p) A A (a|f , <r m , 6, q m t^a/b; q, t;p) \m" - XJ [q -z mt * n -* b / aM . q , t . 



(4.24) 



(4^25) 



If b = q k oi b — t k there is a further symmetry: 
Corollary 4.8. If t(X) < n, 

'k n + X\ C° kn (q~ k ;q,t;p) A^q k a\t n ,t 1 - n q k a;q,t;p) U 



M / [a, q -k];q,t; P ~ ^ (pqaq k ; q, t; p) A\(q 2k a\t n , t 1 ~ n q k a; q, t;p) \x/ [ q k aq -ky qt . p (J " 2(,) 
IfX 1 <m, 

lm k ■ X\ C^ k (t k -q,t;p) A^t- k a\q~ m , q^^a; q, t;p) / fi 



» /[a.t-U-Mv G^.( M a^;g,i;p) A A (^ a | g -™,^i-fc a;g ,i;p)\A/ [t _ fca tfe] ^ t; ; 
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Proof. If we express the binomial coefficient on the left of the first equation in terms of interpolation theta 
functions, we find that complementing /i gives an interpolation theta function that is still evaluated at a 
partition, namely m" — A. Re-expressing the result as a binomial coefficient and applying complementation 
symmetry gives the desired result. The second equation follows by duality. □ 

Since the bulk difference equation (|4.IIJ) can be viewed as a multivariate Jackson summation indexed over 
a skew Young diagram, we also expect there to be a corresponding analogue of the Bailey transformation. 

Theorem 4.9. The sum 

A x (a\b,apg/bf;q,t;p) y A°(a/b\c/b,f,g;q,t;p) /X\ //A 
A° K (a/c\b/c,apq/bd;q,t;p) K ^_ x A°(a/b\l/b,d,e;q,t;p) [0i6] ., itjp W [a/b,c/b];q,t;p 

is symmetric in b and b 1 , where bb'de = capq, bb' fg = apq. 
Proof. Summing 

A"(a|&', apq/b'f; q,t;p) \ - O [a,b>];q,t;p(fj) [ a /b',b/b'}\q,t;p^ [a/b,c/b];q,t;p 

A° (a/c\b/c, apq/bd; q, t; p) ^ A° {a/V\l/V, apq/b, apq/b'f, apq/b'g; q, t; p)A° (a/b\b'/b, pqa/c, d, e; q, t; p) 

(4.29) 

over v gives l|4.28fl : summing over \i gives the same sum with b replaced by 6'. □ 
Remark 1. The bulk difference equation (14. lip is the special case c — b above. 

Remark 2. The rectangular case A = to™, /i = gives Warnaar's conjectured multivariate elliptic Bailey 
transformation, Conjecture 6.1 of as a corollary. One can also obtain a transformation of Schlosser-type 
sums by setting t = q, b = q~ Nl , b' = q~ N2 , and n = 0. 

This induces a symmetry under the Weyl group D^. 
Theorem 4.10. Define a function £l\/ K (a,b:vo,Vi,V2,V3;q,t;p) by 



tl\/ K (a, b:v , Vx,v 2 ,v 3 ; q, t;p) (4.30) 
\p c \/i J ,(p ( l a '/vo,pqa/vi,pqa/v2 > pqa/v3;q,t;p) /\\ / ^\ 

r ^ Ck C° /K (v /c 1 v 1 /c,v 2 /c 1 v 3 /c;q,t;p)- 1 \vl [vqa 2 :Pqac] . q , t . p \K/ [a/c , b/c];q , t ;p 



where c = y/ 'bvoViv 2 v 3 / 'apq. Then Q\/ K is invariant under the natural action of on («o, Ui, v%, V 3 ); in 
particular, 

tt\/ K (a, b:v , Vx,v 2 , v 3 ;q, t;p) = il\/ K (a, b:v , v x , l/v 2 , l/v 3 ; q, t;p). (4.31) 

Proof. Indeed, the symmetry 

ri\/ K (a, b:v Q ,v 1 ,v 2 ,v 3 ; q, t;p) = Cl\/ K (a, b:v Q , ui, l/v 2 , l/v 3 ; q, t;p) (4.32) 

is simply Theorem 14.91 up to a change of variables; since £l\/ K has a manifest S4 symmetry, the D4 symmetry 

follows immediately. □ 
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From the center of D4, we obtain the following transformation, generalizing the commutation relation of 
Corollary 13. 41 ( which corresponds to the case b = b' = 1/q when evaluated at a partition). 

Corollary 4.11. If bed = b' c! d! , then 

2.^ ao (niu mi mt Mh r - a- n y±st>) 



AO/Aim ma. Em- a i- v ) Minium 222. mt \lh c d-a t-v) 

is invariant under (6, c, d, c', c?') 1— > (&', c', d', 6, c, d). 

The case b — q^ 1 gives rise to the following difference equation: 
Theorem 4.12. Define a difference operator D q ^ < ' n \uo:ui:u2,U3,U4 : ;t;p) by 

(D+ (n) (uq-.ui.U2, u 3 , "4; t,p)f)(. ..Zi...) 

TT ejpqtn-^/uolp) 

11 n2<r<5 M""^i;rt ' 



KKn 



n a + flj n n s«<.. .,■**...), 

,<.<„ ,<,<„ e(P9«,/»o;p)«(2?;p) l J£ & 6( Wj ;p) 

where U5 = p 2 q/t n ~ 1 uoUiU2U3U4, and R Zi acts on functions by replacing Zi by its reciprocal. Then 

D+( n \u :uvM2,u 3 ,U4^p)Rl {n \;qV% 1 ,q- 1 /W,q,t]p) = J2 c Xtl R*^(;u 1: uo;q,t-p) (4.35) 

ACmCA+1" 

where 

A ^ ~ pq pq pq pq ,, x \ x / yi.OV) 

l±\(qt Ui/U |t , Mo „ 2 , „ oU3 , „ oU4 , UQU5 ,q,t,p) \A/ [t«-i Ul /„ 0) i/g];q ]t ; p 

One can obtain other identities from our main identities (|4.11(l and Theorem Ij4.9|l by first applying some 
combination of duality and complementation symmetry, specializing the result so that the sum terminates 
regardless of A and/or k, then analytically continuing. 

For instance: 

Theorem 4.13. The interpolation theta functions satisfy the following connection coefficient identity: 

[P^H; a, b'- q, t; p)]P? m ' n) (; a, b; q, t;p) (4.37) 
_ C° mn (t"- l ab,bq,t;p) A,(^\±t n ,q- m ,-£y-,q,t-,p) 



(tn-iaV, b -; q, t;p) A A (^| * , g— , g, i;p) [^, f]; q, t; 



Proof. If we set e = g~ m in the bulk difference equation H4.11JI . the resulting sum will vanish for fi <£_ m n , 
regardless of A, and we may thus replace the range of summation with C m n . If we express the binomial 
coefficients with upper index A in terms of interpolation theta functions, we find that the remaining factors 
involving A are the same on both sides; dividing by those factors and analytically continuing gives the above 
result. □ 
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Applying Corollary 13. 31 gives the following. 
Corollary 4.14. We have the following connection coefficient identity: 

[P;^(;a',b; q ,t;p)}P* x ^ n \;a,b;q,t;p) = ( X ) (4.38) 

W [*^,£](f-W);g,t;p 

Similarly, 

[R^\ ] a',b ] q,t ] p)]R< n \ ] aMqAp) = O (4.39) 

Combining the two identities for interpolation theta functions gives the following result. 
Theorem 4.15. General connection coefficients for interpolation theta functions are given by the sum 

Remark. These connection coefficients generalize Rosengren's construction of elliptic 6-j symbols in |17j. 

Using the connection coefficient formula (|4.39|) , we can extend the special branching rule 13.43fl to a general 
branching rule. 

Theorem 4.16. Interpolation functions satisfy the following "bulk" branching rule. 

R* x (n+k) (...c cel.... vt^ 1 ; a, b; q, i;p) = £ c XK R< n \. . . Zi . . . ;a,b;q,t;p), (4.40) 

K 

where 

ca« = ( A \ . (4.41) 

Proof. Expand the left-hand side in interpolation functions with parameters (v,b), rewrite those as n- variable 
interpolation functions, and convert the parameters back. This gives the coefficients of the branching rule as a 
sum over partitions, which can be summed via the bulk difference equation l|4.11|l . □ 

Similarly: 

Theorem 4.17. The interpolation functions satisfy the following generalized Pieri identity: 

II J {V t!±i P)7 \ R * n) (' • • Zi ■ ■ ■ 5 a > q ~ mh ' q > *5 = E C ^ R « n) (' ■■Zi---;*,b;q,t;p), (4.42) 



wh 



ere 



AO un-l lhU n-l i , v ^{t^alb^^pq/vb^q^/b-q^p) / k\ 
c KX = A mn (t v/b\t va,v/a;q,t;p) ., m ,, u . — — —77- — r ( J • (4-43) 



A A ( g ™t»-ia/6|i», M / u 6, g ™ v /6;g,i;p)\A/ [t „_ la/fc(? _ m] . git . p 
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If we expand a Cauchy-type interpolation theta function in other interpolation theta functions, the result 
can be analytically continued again; we obtain the following result. 



Theorem 4.18. Interpolation theta functions satisfy the following "Cauchy" identity 

A M (ffi|f",g- m ;g,f;p) 
Cl n (t^ab,b/a;q,t;p) " 



A (— -\t n q~ m -q tp) 



l<i<n,l<j<?n fiCm 

K(n,m) 



K^-Ji- ■ ■ Vp/vi ■ ■ ■ ; Vp/«. V*, i/?;p) 
5 Biorthogonal functions 

Definition 12. Let to, ti> ^2, is, web Wij 9, i be parameters such that t 2n ~ 2 totit 2 t 3 uoUi = pq. Then define 

Rx (]t :t 1 ,t 2 ,t 3 ;u ,u 1 ;q,t;p) := V ■— r^rr; - tQTtl ; ; -• (5.1) 

Remark. This is an analogue of the "binomial formula" of |12| . 

The specific form of the above expansion is effectively determined by the requirement that be symmet- 
rical in to through £3; more precisely: 

~ (n) 

Theorem 5.1. The function R\ satisfies the symmetry 

B(n) f , . . , . s R < ^ ) (;t a :t 1 ,t 2 ,t 3 ;u ,u 1 ;q,t;p) 

R\ '(■ 1 ti:t ,t 2 ,t 3 ;uo,ui 1 q,t- 1 p) = — ^y- : -, (5.2) 

R\ '(. ..t n H x .. .■,t :t 1 ,t 2 ,t 3 ;u ,u 1 ;q,t;p) 

where 

R,[ n) (. . . t n -% t :h,t 2 ,t 3 ; u , Ui',q, t; p) = A^l/wo^ilf 1 " 1 ^, t n_1 W 3 , t 1 - n /tiu 1 ,pqt n - 1 t /u ; q, t; p). 

(5.3) 

Proof. If we expand 



R[ n \; t :tx,t2, t 3 ; u , u%;q, t;p) (5.4) 



in terms of R*^ n \; t±, uq; q, t;p) using the connection coefficient formula (|4.39|) . the result simplifies via the bulk 
difference equation l|4.11|l to give the desired identity. □ 

Remark. Alternatively one can apply the same steps used above to derive the connection coefficient identity, 
but starting with Theorem 14.91 instead of l|4.11|l . The resulting computation is somewhat more complicated, 
but does have the merit of deriving the expansion above. 

In particular, we have: 



jW,., , . . ... 1 R*x n \;h,u ;qA;p) 



R x (; t : tl ,t 2 ,t 3 ; u , ^->q, t;p) - ^-^-^ 



(5.5) 



In addition, from the difference equation (|3.34() for interpolation abelian functions, we directly obtain the 

~ (n) 

following difference equation for R x . 
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Lemma 5.2. The function R\ ' satisfies the difference equation 

£><"> (n , t , h, t^p/uotoh; q, t; p)R { ^(; q l / 2 t Q :q l/2 h,q' 1/2 t2, q~ 1/2 t 3 ; q 1/2 u , <T V V; <?, t; p) (5.6) 

= Y[ O(t n - l u t , t n - l u t u e-Hohip)^ (; Mi, t 2 , is; "o, tti; g, t; p) 

l<i<n 

Theorem 5.3. The functions R^ agree with the biorthogonal functions of 115$ . To be precise, in the notation 
of that paper, 

to-ti,t 2 ,t 3 ; u , ui;t;p, q) = R^/{; Mi, h, h; u , m;p, t; q)R^(; Mi, t 2 , *3; u , m; g, (5.7) 
Proof. We need to show 

^o"a(; Ml, *2, *3! uq, Ui;t;p, q) = R^" () Mi, *2) *3! "o, «i! <7, *; P)- (5-8) 
This equation is preserved by the action of the difference operators 

D (n) (uo,ir,i s ,i 1_ >/u M s ;g,i;p), < r < s < 3; (5.9) 
in particular, both sides satisfy the same generalized eigenvalue equation with respect to the compositions 

D<*>(uo, to, ti;q, f;p)£< n V /2 «o> <T 1/2 * 2 , <T 1/2 i 3 ; 9, t;p), (5.10) 

D(")( M o,to,t 2 ;g,i;p) J D ( " ) (? 1/2 «o,^ 1/2 ti,^ 1/2 i3;g,t;p),... (5.11) 

and thus agree up to scalar multiples. Since both sides evaluate to 1 at . . . tot™~ 1 ■ ■ ■ , they agree everywhere. □ 

Remark. Similarly, one can show using Theorem l4.12l that the difference operator defined there acts as a raising 
operator on the biorthogonal functions, just as in |15|. 

In |15| , it was shown that these functions are biorthogonal with respect to an appropriate contour integral, 
and a number of other properties were given. A few remaining properties were outside the scope of that 
paper; we are now in a position to prove these. Most striking of these is evaluation symmetry, generalizing 
the analogous property of Macdonald and Koornwinder 03 1^| polynomials (see also ^B] for a proof for 
Koornwinder polynomials along the present lines). 

Theorem 5.4. For otherwise generic parameters satisfying t 2n ~ 2 totit2t 3 uoUi = pq, 

^(...tof-V ...;t :t 1 ,t2,t 3 ;u ,u 1 ;q,t;p) = R[ n \. . .tot^q^ . . .;i :ii,i 2 ,i 3 ;uo,ui;q,t;p), (5.12) 
where 

I ; ~ - - - ~ „ Uq Uo Ui Ul 

to = \Jt tit 2 t 3 /pq t ti — t a ti t t 2 =t a t 2 t t 3 = t t 3 — = — — = —. (5.13) 

to to to to 

Proof. Upon specializing the variables as required, the result can be expressed as a sum over binomial coef- 
ficients; changing from the original parameters to the "hatted" parameters interchanges the binomial coeffi- 
cients. □ 
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Remark 1. Aside from some simple factors, this sum over binomial coefficients is the analytic continuation of 
the "multivariate 6-j symbol" of Theorem 14. 151 

Remark 2. When the biorthogonal function is specialized to an interpolation function, we obtain the identity 
R; {n) (. . . (u/a)i"- V* b/a; q, t; p) #* (n) (. . . {v/a'^q^ . . . ; a', b/a'; q, t; p) 



#* (n) (. . . (w/o)*"-* . . . ; a, b/a; q, t;p) R* x {n) (. . . {v/a')t^ . . . ; a', b/a'; q, t;p) 

where 



(5.14) 



a' = v / t n ~ 1 bv/a. (5.15) 

If v = q~~ m t x ~ n here, we recover the complementation symmetry of binomial coefficients; by analytic continua- 
tion, the two results are equivalent. 

Remark 3. In the univariate case, this was proved in Section 9 of |27| . 

Using inversion of binomial coefficients, one can expand interpolation functions in terms of biorthogonal 
functions. 

Theorem 5.5. Biorthogonal functions satisfy the following "inverse binomial formula". 

tliM/ 4-4-4-4- 4 MD*(n)/ 4 4 \ W [t™ 1 1 / ' U ,t™ 1 1 U 1 ] j <J , t \p „N 

[Ryhto:ti,t2,t a iUo,ui\q,t\p)]R^ ' (;t ,u Q ;g,t;p) = - 1{ — pq pq — — (5.16) 

Combined with the binomial formula (|5.1|) . we obtain connection coefficient formulas, analogous to connec- 
tion coefficients for Askey- Wilson polynomials 

Theorem 5.6. If t 2n ~ 2 tQtit2tzUQUi = pq and t'^t'^u^ = tit2t^ui, then 

[^1™K;*o:*i,4>^; w oXi;3,*;p)]-R> (\ta-ti,t2,t3;uo,ux;q,t;p) (5.17) 



~ K( ^-i x A° (i^ajtn-itot!, t^-itota, ^-itot 3 , t"-H oUl ;q, t;p) 



^ [^H' tn-i t0M1 ];g. f ;P ^ [t™- 1 *o/«o,*™- 1 to«i];g,t;p 

If t' 3 = t 3 , the same connection coefficients can be computed via R*( n >(; t 3 , uo; q, t;p); the result is precisely 
our generalized Bailey transformation, Theorem l4.9l If also t' 2 = t 2 , the bulk difference equation (|4.11() applies, 
giving the following result announced in |15| . 

Corollary 5.7. 

[RlT ) (;to-tiV,t2,t 3 ; u 0l ui/v;q,t;p)}R[ n \;t :t 1 ,t 2 ,t 3 ;u , Ux;q,t;p) (5.18) 

\/V [l/uoui,l/i)](t"- 1 t 2 t3,piit™~ 1 to/«o,ti'"/«i);'3,t;p 

We also have discrete biorthogonality, which was derived in |15j via residue calculus, but can also be derived 
via the present "hypergeometric" methods. 
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Theorem 5.8. For any partitions A, k C m n , and for otherwise generic parameters satisfying t^ti = q m t n , 

t n -H 2 t 3 U0Ui =pq m+1 

Rx ] (• ■ ■ W'V' t :t u t 2 ,t 3 ; u , m; q, t;p)R<f (. . . t Q t n -' l q^ . . . ; t :h,t 2 , t 3 ; Ul , u Q ; q, t;p) (5.19) 

A^(t 2n ~ 2 t 2 \t n , rVi.t"" 1 ^, t n -H t 3) t n -H u , t n -H u 1]q , t;p) = 

unless X — k, when the sum is 

ao ( f-Hx iti _pq_ _e_ _pq_. n f . n \ 



(5.20) 



Proof. The argument of |16j carries over essentially verbatim; alternatively, the argument of |17| generalizes, 
using Theorem 14 .151 □ 

Remark. Note that the above inner product is normalized by equation (|4.15|) . 

We also have a special quasi-branching rule (having the branching rule for interpolation functions (Theorem 
14.161) as a special case): 

Theorem 5.9. The biorthogonal functions satisfy the expansion 

R^ l+k \. . .Zi. . . , t ,t t, . . .t Q t k ~ 1 ]t :ti,t 2 , t 3 ;u , ui;q,t]p) (5-21) 

= ^ c ^ R< k ] (■■■Zi---; t t k :ti,t 2 ,t 3 ] u , uit k ; q, t;p) 

where 

c Xk = ( X ) +k i (5.22) 



Proof. Expand the left-hand side in interpolation functions, apply equation (|3.43[> . and expand back into 
biorthogonal functions; again the resulting sum can be simplified via the bulk difference equation (|4.11|l . □ 

Similarly, we obtain a special quasi-Pieri identity. 

Theorem 5.10. The biorthogonal functions satisfy the expansion 



9(t Q z l ,t /z i ;q;p) m - („) 

• ■ • ■ ■ ■ ,q 

^ c KA ^l n) (- - -Zi-. .;t Q :ti,t 2 ,t 3 ;u ,ui;q,t;p), 



— — - — — - — —R\'{...z 1 ...;q t :ti,t 2 ,t 3 ;q u ,Ui;q,t;p) (5.23) 

i<i<n e \\PQ/ u o) z i> {pq/ u o)/Zi]q;p)m 



where 



t n ~ ^ t 

CkA = A° mn ( ^t n -H tx,t n -H t 2 ,t n -Hot 3 ,t n -Houi]q,t]p) 

A / [_i_,g-™](tn, 1 ^a);g, t;p 

- u o u i n ' t™- 1 *oui ' u o ' r 
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Finally, we have a Cauchy identity for biorthogonal functions. 
Theorem 5.11. The function 

Fix, x n ; Vl y m ) := ni<,<„,K 3 <Jfe^^ 3 M;p) 

ni< 4 <n% m u Xi,q m u /x i ]q;p) m ]l 1 < j < m 9(pq m /u y j ,q m y j /uo]l/t;p) n 

admits an expansion 

t n u i n_ V 



c ^ n) (- -- x i-- ■;t a :t 1 ,t 2 ,t 3 ;uo,u 1 ;q,t;p)R ( l ™ ) ^ fl ,(. ..y t .. . ; t :ti, t 2 , t 3 ; -—± ;t,q;p), 

C^ nn (pq/uQto,q m t /u -,q,t;p)F(xi, . ..x n ;yi, . ..y m ) 
~ A° mn (^-Ho/uol^-Hah^^Hoh^^Hots, t n -H QUl ;q, t;p) 



(5.26) 



where 



1 1 n m t 

c^ = A^ \t n ,q- m ,g_ £ ;g>f;p ). (5.27) 

U Ui t n HqUi u a 



6 Algebraic modularity and rationality 

The purpose of the present section is to give a purely algebraic definition of the biorthogonal functions; as 
a consequence, it will follow that the biorthogonal abclian functions are in addition modular functions. This 
could of course be shown directly by determining how the interpolation functions and C* symbols behave under 
modular transformation; the approach used here has the advantage of being more conceptual in nature. In 
addition, we also obtain analogues of interpolation functions, binomial coefficients, and biorthogonal functions 
over arbitrary fields (including those of positive characteristic). 

In addition to making sure the construction is sensible geometrically, wc also want things to be reasonable 
arithmetically; that is, in such a way that the functions depend rationally on the parameters. In particular, 
the difference equation as given above is problematical in this respect, as it requires us to choose a square root 
of q, despite the fact that the interpolation functions themselves are independent of that choice. The simplest 
way to fix that is to relax the notion of elliptic curve slightly, by forgetting which point represents the identity; 
this allows us to absorb the freedom in choosing the square root. We thus obtain the following fundamental 
definition. Recall that for a (smooth) curve C, the Picard group Pic(C) is the group of divisors modulo principal 
divisors; for an integer n, Pic n (C) is then the preimage of n under the natural degree map. For each n, Pic™(C) 
has a natural structure of algebraic variety, a principal homogeneous space over the algebraic group Pic°(C). 
We will use multiplicative notation for divisors (and the group law in Pic(C)), and will denote the divisor 
associated to a point p € C by (p), and to a function / on C by (/). 

Definition 13. A genus 1 hyperelliptic curve is a pair (C, r), where C is a genus 1 curve, and r e Pic 2 (C) 
is a divisor class of degree 2 (which thus induces an involutory automorphism t T : x i— » t/x of C with genus 
quotient). A BC„(r)-symmetric function on C is a function / on C n invariant under permutations of the 
variables and replacements Xi \— > r/xi. 

Remark 1. Note in particular that Pic°(C) is an elliptic curve, and Pic^C) can be canonically identified with 
C. Each point x € C thus induces a map Pic(C) — > Pic°(C) defined by y i— » yx~ des( - y \ and in particular gives 
an identification of C with the elliptic curve Pic°(C). 
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Remark 2. The action of l t on C extends to an automorphism of the Picard group Pic(C), which we will denote 
by L r (x) :— T dcg ( x ) jx. Similarly, its action on the zth copy of C in C™ will be denoted by n^ T . 

Remark 3. For x € Pic°(C), there is a natural map from the genus 1 hyperelliptic curve (C, r) to the curve 
(C,x 2 t) given by p ^ xp. Since all of the functions we will define below are canonically defined, they in 
particular will transform nicely under this isomorphism; in general, replacing each element y e Pic(r) by 
x de s(v)y w iH leave the function invariant. 

While there is in general no canonical choice of divisor representing r, there is in fact a canonical choice of 
representative for r 2 , namely the ramification divisor of the quotient map C i— > C/i T . In odd characteristic, this 
is simply the product of the fixed points of i T in Pic (C); in even characteristic, it is the square or fourth power 
of that product (depending on whether the curve is ordinary or supersingular) . We will denote this divisor by 



Given a divisor D, let C{D) denote the space of functions / on C such that (f)D is effective; by the 
Riemann-Roch theorem, this space has dimension max(deg(_D), 0), unless deg(D) = and D is principal 
(when the dimension is 1). The following proposition is key to proving the difference equation in arbitrary 
characteristic. 

Proposition 6.1. Let D be a divisor of degree m. Then l* t C(Dl t {D)) = C(Dl t {D)), and the dimension of the 
space of t* -invariants is max(m + 1,0). Moreover, 



with images of dimension max(m — 1,0), max(ra + 1,0), respectively. 

Proof. The statement about the space of t*-invariants follows immediately from Riemann-Roch, since i*- 
invariants are simply functions on the quotient P 1 . For the remaining statements, multiplication by a l*- 
invariant function allows us to replace D by any other divisor of degree m; in particular, we can assume that 
D is supported away from the support of (t 2 ). That 



then follows from the definition of ramification divisor; the dimension of the image can be computed from the 
dimension of the kernel. The remaining claims follow upon multiplication by a t*-anti-invariant function. □ 

Fix a genus 1 hyperelliptic curve, as well as two points q,t £ Pic°(C). Given four generic points uq, u\, 
u 2 , u 3 e C such that t n ~ 1 UoUiU2U 3 = r 2 , we define a difference operator D^ n \;uo,ui,U2]q,t;C,T) acting on 
-BC„(r/g)-symmetric functions as follows: 



(D^ (; n , ui.ua; q, t; C, r)/)(. .. Xl ...)= \\ (1 + < T )<5 (n) (- • • a* • • • ; uo, «i, «2, «s; *; C, r)/(. . . Xi . . . ) (6.4) 




(l-t;)£(D, T (D))c£(Dc T (D)/{r 2 )), 
(1 + l*)£(Dl t (D)(t 2 )) c £(Dl t (D)), 



(6.1) 
(6.2) 



(1- l* t )C{Di t (D)) c C{Dl t {D)/{t 2 )) 



(6.3) 



l<i<n 



for a function 5^ defined inductively as follows: 
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(0) 5°{u a ,u 1 ,U2,u 3 ;t;C 1 T) = 1. 

(1) As a function of x n , S^(. . . Xi . . . ; uq, u\, v,2, u 3 ;t; C, r) has divisor 

l<i<n ( tXi ) 



(6.5) 



(r 2 )ni<i<„<^) 

(2) Setting x n = uo gives 

5 (ti) (xi, . . . ,x„_i,u ;wo,wi,M2,u 3 ;t) = 5 (n_1) (xi, . . .,a; n _i;iiio,«i,«2i«3;t) (6-6) 

Proposition 6.2. For fixed (C, r) andt, 5^ n ' is a well-defined rational function of all variables and parameters. 
It is invariant under permutations of X\ through x n and under permutations of Uq through u 3 . 

Proof. For n — 1, we find from Proposition 16 . II that 

(1 + L*)S ( - 1) (x;uo,ui,u 2 ,u 3 ;t;C,T) (6.7) 

is a constant; evaluating this constant at x = u$ shows that it is 1. Evaluating at x = u% shows that 

S ( - 1 \x;u ,ui,u 2 ,u 3 ;t;C,T) = S (1) (x;ui,u ,u 2 ,u 3 ;t;C,T) (6.8) 

as required, and thus the proposition holds when n — 1. 

For n > 1, we first observe by induction that 6^ n ' is symmetric in x\ through x n -\ and in ui, u 2l u 3 . Now, 
if we set x n — uq, the divisor in x n -i will be 

\ tuo 



(r 2 ) Ill<i<n-l(xi/ 



(6.9) 



from the dependence on a;„_i of the divisor in x n , we deduce the dependence on x n of the divisor in x n -%, and 
thus conclude that the divisor in x n -\ of 5^- n ' is 

{^i(^i(ti2i(u3> rii#n— 1 \ txi i ir, -. ,-,\ 

We claim that setting x„_i = u\ gives , u 2 ,u 3 ) = <5 ( ™ 1) (;tMi,u ,U2,U3); indeed, setting x„ = Mo, 

x n -i = u\ gives S^ n ~ 2 \;tuo,tui,U2,u 3 ), and the divisor in x n is also correct. It follows that 

S (n) (x 1 , . . .,x n ;u ,ui,u 2 ,u 3 ;t;C,T) = S (n \x 1 , . . . ,x n -2,x n ,x n -i; ui,u , u 2 ,u 3 ; t; C, r), (6.11) 

at which point the proposition follows. □ 

Remark. If C is a complex elliptic curve of the form C*/(p>, then can be expressed in theta functions as 

5 (n) (. ..Xi...;u , ui,u 2 ,u 3 ;q, t; C, r) (6-12) 

TT 9(u Xj/T, UlXi/T,U 2 Xi/T, U 3 Xj/T;p) jj 0(tXiX 3 /T;p) 

~ J t { n 9{x 2 jT)e(u 0Ul t n -yT, U U 2 t n -y T , Ua U 3 t^/T-p) x< l} j<n 0(x iXj /T; P ) ' 

where the parameters are lifted to C* so that uqU\U2U 3 = pr 2 ; thus the difference operator defined above agrees 
(up to a scale factor) with the earlier definition. 
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For an integer m > 0, let A m (u ; q; C, r) C ^■(Y[i<i< n ( < l^ u o) ( q -i Uo )) denote the subspace of t r -invariant 
functions, and let Am (uq; q; C, r) be the n-th symmetric power of that space, viewed as a space of functions 
on C"; in particular, Am) [uq] q; C, r) consists of BC n (r)-symmetric functions. We also write A^(uo;q;C, r) 
for the union linim^oo im' (no; 9; C, r). 

Theorem 6.3. For any integer m > 0, 

£( n ) («o, ui, u 2 ; q, t; C, r)A^ (u ; q; C, r/q) C i4&> (u ; g; C, r). (6.13) 

Moreover, 

{D in+k) (u ,u 1 ,u 2 ;q,t;C,T)f)(z 1 , ...,z n ,u ,.. . ,t fc_1 -uo) 

= (Z)< n > (u , «i, «25 9, *5 C, r)/(_, u , . . . , i* -1 ^))^ . . . , z n ). (6.14) 

Proof. The second claim follows readily from the definition of the difference operator and of 8^ n > , so we need 
only consider the first claim. In particular, it suffices to show that as a function of z n , the only poles of 

^ (ri) ( M o,wi,M 2 ;(?,i;C,r)/ (6.15) 

for / £ Am\u a ; q; C, r/q) are of the form q~ l u$ or qV/u f° r 1 < i < m, of multiplicity at most 1. Now, as a 
function of z n , 

6 (ti) (. . . Zi . . . ; u , ui,u 2 ; q, t; C, r)/(. ..z t ...) (6.16) 

has polar divisor at most 

(r 2 ) II II <<r^><<r 1 r/«o>/<T/uo>, (6-17) 

l<z<n l<i<m 

with the last factor coming from the numerator of S^ n \ In particular, upon symmetrization by t T , the result 
clearly has polar divisor at most 

(r 2 ) H {zi){T/Zi) J] (9 _, «o)<?V«o)- (6-18) 

l<z<n l<i<m 

We thus need only show that the potential poles corresponding to the first two sets of factors disappear upon 
symmetrization; we may assume (since it is true generically) that the three sets of factors have disjoint support. 

That the singularities corresponding to (r 2 ) disappear follows immediately from Proposition ^. II so it remains 
by symmetry to consider the potential singularity at z n — z n _\. Now, in the sum 

Y[ (l + il T )S^\...z l ...;u ,u 1 ,u 2 ; q ,t;C,r)f(...z l ...) (6.19) 

l<i<n 

= I[ ( 1 + L lr)[( l + 4-l,r + Cr + Ci.tC)^"' (...Zi...;U , U X ,U 2 ] q, t) C, t)/(. . . Z t . . . )] 

\<i<n-2 

the factors corresponding to 1 and i* n _i T i* n T are already nonsingular at z n — z n -i, so it suffices to show that 

(Ci,r + C,rM (n) (...Zi...;uo, Ul ,u 2 ; q, t; C, r)/(. . . z t . . .) (6.20) 
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has no pole at z n -\ = z n . Now, since 5^ and / are symmetric under permutations of z\, . . . , z n , this sum can 
be written as the symmetrization of 

£ iT <f<"> (. . . Zi . . . ; u , u u u 2 ; q, t; C, r)/(. . . z { . . . ) (6.21) 

under exchanging z„_i and z n ; the disappearance of the pole then follows from local considerations. □ 

We will also need to construct algebraic analogues of certain products of C* -symbols. First, we need the 
following family of univariate functions. For a nonnegative integer to, let u , . . . , u 2m -i eCbea collection of 
2to points such that u ui ■ ■ ■ u 2m _i = r m , and let I? be a divisor of degree to. Then there exists up to scale a 
unique function / with divisor 

n <i<2 TO <«i> (6 22) 



(Dir(D)) ■ 

Thus the function 

fl x ) 

uj(x:u , . . . , u 2m -i;C, t) := (6.23) 

f(r/x) 

is well-defined; multiplication of / by an i*-invariant function shows that uj is independent of the choice of D. 
Note that 



Lemma 6.4. The function w has divisor 

llo<i<2m( U «) 



(6.24) 



and satisfies the normalization that for any point xq G C such that x\ = t, 

uj{x q :uq, . . . ,u 2m -i;C,T) = 1. (6.25) 

Moreover, 

u(x:u 0l . . . , u 2m -i, v, t/v; C, t) = u>(x:u 0l . . . , u 2m -i; C, t), (6.26) 

and 

lo(x:u , u 2m -i.C, t)uj(x:v , v 2n -i;C, r) = lo(x:u , u 2m -i,v , v 2n ^C, r). (6.27) 
Proof. The first claim is straightforward; for the second claim, we find that 

oj{x:u ai . . . , n 2ro _ i; C, r) - 1 = f{x) ' ~ f . { \ ,x) G ^((r 2 )- 1 ), (6.28) 

and thus it vanishes as required. The final claims follow immediately from the original definition of w. □ 

Another univariate function can be obtained from the observation that since r induces a map to P 1 , any 
four points induce a uniquely defined scalar, namely the cross-ratio of their images. We denote this by 

xia^^J^Cr), (6.29) 

which can also be defined as the unique i*-invariant function of a such that 

X{b, b, c, d; C, r) = 0, x(c, b, c, d; C, r) = 1, x(d, b, c, d; C, r) = oo. (6.30) 
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We can also define this by 

X(a, b, c, d; C, r) = ui(x:bax/T, bx/a, dx/c, dcx/r; C, x 2 bd/r), (6.31) 

for any x £ C. Note also that since \ IS defined as a cross-ratio, all of the usual transformations of cross ratios 
apply; note in particular the symmetry 

X(a, b, c, d; C, r)%(a', b, d, d; C, r) = x(a, b, d, d; C, r)x(a',b, c, d; C, r). (6.32) 

In addition, the identification with u>, with its different natural symmetry, gives rise to the transformations 

X(a, b, c, d; C, r) = x( a , b, c, d; C, abed/r) (6.33) 
ui(x:v ,v 1 ,v 2 ,v 3 ;C,t) = uj(x:v q , vi, x 2 /v 2 , x 2 /v 3 ; C, tx 2 /v 2 v 3 ). (6.34) 

Definition 14. Let a, vo, ■ ■ ■ , v 2m -\ € Pic°(C) be such that 

v vi---v 2m -i =q m a m . (6.35) 

Then for any partition A, we define 

A° x (a\v , ■ ■ ■ ,v 2m -i;q,t;C) = [[ uj(q 1 ~H' 1 ^xwqx, . . . , v 2m -ix; C, qax 2 ), (6.36) 

for a generic point x e C (on which the value docs not depend). 

Definition 15. Let a, bo, . . . , b 2 i-i € Pic°(C) be such that Ilo<r<2i b r = (t/q)(qa) 1 ^ 1 ■ Then we define a symbol 
A\(a\b 0} . . . , b 2 i-i;q, t; C) for all partitions A via the recurrence 



■f n+1 c 

A x (t- 2 a\b a /t, . . . , bu-i/t, q^/t, q~ m , q m+1 a/t, q m a; q, t; C) 
with 

A (a|6o,-.., 621-1; g,*;C) = 1; (6.38) 

here, of course, the limit is to be taken in the algebraic sense: the evaluation at the appropriate point of the 
stated rational function in the limit variable. Similarly, we define ( x )^ a b y q t . c via the recurrence 

™'V) = lim A° n (a/b\q m a/x,l/bx,qax/b,qx/q m ;q,t;C) (6.39) 

A° x {a/t 2 \b, qa/tb, q^/t 2 , q m a/t; q, t; C) (X 



A m .x(a\b , . . . ,b 2 i-i;q,t;C) = lim A m (a|fe , • • • , foi-i, q m+ a/x,qa/xt,qax,qx/q m ;q,t;C) (6.37) 

x— »1 



A° x (a/bt 2 \l/b, qa/t, q^/t 2 , q™a/tb; q, t; C) \\J [a/t , M , q , t ,c 

with 

= 1. (6.40) 



[a,b];q,t;C 
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Remark. For a complex elliptic curve, the symbols A° and A\ are precisely those defined above; the parameters 
must be lifted so that the constraints on Yii bi are satisfied as given, except with q replaced by pq. The remaining 
symbol corresponds to the product 



C x (a;q,t;p)C x ) (l/b,pqa/b;q,t;p) ( A 



C+ {a/b; q, t; p)C° (b, pqa; q, t; p) \\J [a . b] , q ^, p 
The algebraic analogue of interpolation functions are defined as follows. 



(6.41) 



Definition 16. For generic parameters a, b, v £ C, q,t £ Pic°(C), we define the algebraic interpolation function 
R* x (; a, b(v);q, t; C, r) as follows. 

(1) R*W(;a,b(vy,q,t;C,T) £ A^(b;q;C,r). 

(2) For each integer m > Ai, and every partition /i C m" with /i ^ A, let I be as in the definition of p*( m - n ) 
above. Then for generic c £ C, 

( II X(zi, q- J b, t n ~*v, c; C, r)< n) (; a, b(v); q, t; C, r)) (bq-^ , . . . , fccT^ 1 , a^f 1 "'" 1 , . . . , aq^) = 0. 

l<?<n 
l<j<m 

(6.42) 

(3) R* x n ^ satisfies the normalization 

R* x {n) (. ..vt n ~ l ...;a,b(v);q,t;G\T) = 1. (6.43) 

As before, these functions are generically uniquely determined by the conditions; moreover, they satisfy a 
difference equation, the proof of which is a direct analogue of the corresponding proof for theta functions. 

Theorem 6.5. The functions R* x satisfy the difference equation 



D {n) (a, b, c; q, t; C, r)i^ (n) (; a, b(v);q, t; C, r/q) (6.44) 

t n ~ 1 a avt 71 ' 1 acqt 71 ' 1 r qr . n s u *{n), u x . n \ 
A A ( — 7— I > ' 7T' T~'^ C ) R X \\a,b{v);q,t;C,T) 

T T OV OC 



and the normalization condition 

R* x {n) (. . . v'f 1 - 1 . . . ; a, b(v);q, t; C, t) = ^{t^a/b^-^av' /t, a/v', qr/vb, f^qv/b; q, t; C). (6.45) 

Proof. As mentioned, the proof follows that of Theorem 13. II the only difference is the treatment of the normal- 
ization. We find as before that 

D^(a, b, c; q, t; C, r)R x {n) (; a, b(v);q, t; C, r/q) cx i?* (n) (; a, b(v); q, t; C, r); (6.46) 

for c = v, we find immediately that 

B W (a, b, v; q, t; C, r)R* x {n) (; a, b(v);q, t; C, r/q) = R* x (n) (; a, b(v);q, t; C, r). (6.47) 
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Now, if we evaluate at the point (. . . at n l q Xi . . . ), we find that only one term survives, and thus the constant 
of proportionality in the difference equation can be computed as 

$W (. . . ai n -*g A< . . . ; a, 6, c, T 2 /abc; t; C, r) 



£(«) (. . . at n - l q x i ...;a,b,v, T 2 /abv; t; C, t) 
We readily verify from divisor conditions that the function 
£(")(. . . Xi . . . ; a, b, c, f-^/abc; t; C, t) _ 

Ki<n 



(6.48) 



= J] x{^lcM n -\rlv-C,t n ~ x ab), (6.49) 



£(")(. . . . . . ; a, 6, v, t x - n T 2 /abv\ t; C, t) 

and we thus obtain 

(5 (n) (. . .at n - l q x ' . . .;a,b,c,T 2 /abc;t;C,r) -p-r , A . ( x 

^)(...a^...;a,M,rVa6,;^,,) = JI*^ 'A^^" \t/v;C, t n x ab) (6.50) 

On the other hand, we find 

A^(r- 1 a/6|a^™- 1 /T,acgt"- 1 /T,T/6w, g r/&c;g,i;C) = J] J] x(ai"~Y, t/c, at^q^ 1 , t/v; C, f^afe) 

l<i<n l<j<Ai 

(6.51) 

= II xK^^.^fl^VMC.r'tti)) (6.52) 

l<i<n 

as required. 

The formula for changing the normalization follows by applying the difference equation to both sides of the 
equation 

Rl {n \ ] aXv)-,q,t ] C,T)=Rl {n \ ] a,b{t 1 - n T/v)-q,t ] C,T), (6.53) 
which gives a recurrence for the relevant scale factor, with the stated solution. □ 

In particular, we obtain the extra vanishing conditions just as before. The various identities of Section |3| 
for the analytic interpolation functions all carry over to the algebraic interpolation functions. In particular, we 
note the Cauchy case 

R* x {n \...z t ...;a,b(v);q,t;C,t n ab/q)= j [ at n_A V _1 , «* n_i , a*V _1 ; C, t n ab/q) (6.54) 



l<i<? 
1<? 



Proposition 6.6. We have the identities 



i? m „_ A (. g, t; C, r) = n — j^^h^T ^ c> r) , (6-55) 

l<j<rn 

C+A • ■ ■ ^ ■ ■ ■ ; a, b(v); q, t; C, t = _ , 6.56) 

1 <j<m 

i?* ( " +fc) (. . . «j . . . , a, ai, . . . , at k ~ 1 ;a, b(v);q, t; C, t) = A° (t n+k - 1 a/b\t n , qa/bt, t k a/v, qt n+k - l v/b\ q, t; C) 

Rl {n) (...z l ...;t k a,b(v);q,t;C,T), (6.57) 
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and thus 

( X ) 

i? A (...at g ..., a ,6( W ),a,t,G,r)- ^^^^/^^n-^^^.^cy ^ 

Proof. The first two identities are straightforward, as both sides satisfy the same vanishing conditions and 
normalization. For the second identity, the vanishing conditions are trivial to verify, and thus both sides 
are proportional. From the difference equation, it follows that the scale factor is independent of r, and we 
may thus reduce to the Cauchy case, for which the verification is straightforward. The final equation follows 
immediately. □ 

We also have the following symmetries which are trivial consequences of the isomorphism invariance of our 
definitions: 

R*^ n \. . . Zi . . . ; a, b(v); q, t; C, r) = R^ n \. . . (x/a)zi . . . ; x, xb/a(xv/a);q, t; C, x 2 t/o 2 ) (6.59) 

= R\ {n \. .. Zi ...; r/a, t/Ht/v); l/q, 1/t; C, r) (6.60) 
= R^ n \. bir/e^v); q, t; C, r) (6.61) 

With the above in mind, we define the algebraic binomial coefficients as follows. For a genus 1 curve C and 
a,b,q,te Pic°(C), we define 

A J := A„(£|t n , 1/6, v, t x - n a/v; q, t; C)R* l } n) {xq x H 1 - l ;t 1 ~ n x, xb/a(xv/a);q, f, C, x 2 /a), (6.62) 

AV [a,b];q,t;C " 

where the right-hand side is independent of n, x e C, v € Pic°(C). (In particular, the binomial coefficients 
really are elliptic in nature.) When /i = A, this is of course consistent with our previous notation Q)^ ^ t _ c . 
If C is a complex elliptic curve, then we obtain the analytic binomial coefficients: 

' A ) =( A ) ■ (6.63) 

AV [a,b];q,t;C/{p) VAV [a,b];q,t;p 

Similarly, if v viv 2 = q 2 a 2 /b, we define 

'A\ A° x (a\b,v ,v 1 ,v 2 ;q,t;C) ( A 



, (6 64) 

{aM{vo , VuV2)w ,t;C &°(a/b\l/b,V ,V U V 2 ;q,t;C) W [a,6];,,t;C 

In particular, the same proof as in the analytic case gives the algebraic bulk difference equation. 
Theorem 6.7. For otherwise generic parameters on Pic°(C) satisfying bcale = aq, 
A\ A£(a/c|l/c,6d,6e,aa/6;<?,t;C) ^ /-n^ A \ /V\ 

W [a,c\; q ,t;C A^C, bd, be, aq/b; q,t; C) " VaV [a,&];g,t;C W [a/6,c/6];g,t;C 

(6.65) 

In particular, 

Q = (6.66) 

„ \AV [a,6];g,t;C W [a/b,l/b];q,t;C 
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Remark. One immediate consequence is Spiridonov's observation |25| that the various more traditionally hy- 
pergeometric special cases discussed following Theorem 14.11 above are modular and abelian in all parameters. 
(This includes the BC n -type sum of ^H|j since as observed there, that sum can be obtained by specializing 
Warnaar's Schlosser-type identity.) 

Theorem 6.8. The algebraic interpolation functions satisfy the connection coefficient identity 

[R;( n H;a\b(v); q ,t;C,T)]Rl (n \:a 1 b(v);q,t ] C,T) = ( X ) . (6.67) 



We define algebraic biorthogonal functions by: 



, r t n , Ut . c K (n) (; *o, «o(*i); q, t; C, t) 

rV>{; t :h,t 2 , t 3 ; u , u V , q, t; C, r) := £ ^ ' t L- ltot3 t ,- HoUl qt ,- Hl _ , (6-68) 

where to, ^2: ^3 , uq, «i € C such that t 2n ~ 2 t tit 2 t 3 uoUi = qr 3 . In particular, when C is complex, we exactly 
recover the analytic biorthogonal functions, which must therefore be modular and abelian in all parameters. 

Theorem 6.9. The algebraic biorthogonal functions satisfy the symmetry identity 

5(n)/ . , , , . n s R ( ™ ) (;t :t 1 ,t 2 ,t 3 - 1 uo,u 1 ;q,t;C,T) 

R\ '(;ti:t ,t2,t 3 ;uo,ui;q,t;C,T) = — ^- : -, (6.69) 

R\ '{. ..t n .;t Q :t 1 ,t 2 ,t 3 ;u ,u 1 ;q,t;C,T) 



wh< 

6(«)/ 4.n-U ++ + + \ aO/ t *1*2 t n ~ 1 tit 3 T qt n ~ 1 t Q 

R\ '{...t t 1 . . .;t :t 1 ,t 2 ,t 3 ;u ,u 1 ;q,t;C,T) — A A ( 1 , , ——r- , ;q,t;C), 

U(jUl T T V t\U\ UQ 



(6.70) 

and the difference equation 

£> (n) (u , t , h;q, t; C, t)R ( ^ (; t :tx,t 2 /q, t 3 /q; u Q , Ui/q; q, t; C, r/q) = R^(; t :tx,t 2 , t 3 ; u , Uf,q, t; C, r) (6.71) 

From isomorphism invariance, we obtain the identities 

jj^ n) (. . . Zi . . . ;t :ti,t 2 ,t 3 ] u ,ui; q,t; C, t) = R^\. . . xz { . . . ; xt :xti, xt 2 , xt 3 ; xu ,xui;q,t; C,x 2 t) (6.72) 

= R^\. ..Zi.. .;T/t :T/t 1 ,T/t2,r/ts;r/uo,T/u 1 ;q,t;C,T). (6.73) 

The statement of evaluation symmetry for the algebraic biorthogonal functions requires a certain amount 
of care, since the usual statement involves a square root; again, this can be absorbed into the hyperelliptic 
structure, but there does not seem to be as natural a way of doing so. What we find is that the hatted 
parameters (including a f ) must be related to the original parameters by 

ti _ qr_ t2 _ qr_ h _ q]_ uo _ «o u\ _ u\ i| _ hhhh , g _ 

i hh ' to hh ' i hh 'to h' to h' f qr 2 

In particular, we can choose to arbitrarily, at which point the remaining hatted parameters are determined. 
This freedom of tg comes from the translation symmetries of C, and corresponds to the fact that the evaluation 
of a biorthogonal function at a partition is really a function on Pic°(C) 7 . 
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The remaining identities satisfied by the analytic interpolation and biorthogonal functions all carry over to 
the algebraic case straightforwardly; for instance, we have the connection coefficient 

[R^(; t :hv, t 2 , t 3 ;u , Ui/v; q, t; C, t)]R^ (; t :h, t 2 , t 3 ; u , u^q, t; C, r) (6.75) 

'A 



\M/ [ T 11/ t" 1 t 2 t 3 1 tQ Vt\ \ f.^, 

'«0»1 '« JV t ' "0 ' "1 ' 

and so forth. 



7 Elliptic bigrids and degenerations 

As we mentioned above, there is a fifth class of perfect bigrids, namely elliptic bigrids. In contrast to the earlier 
classes, the space of elliptic bigrids is of bounded dimension (8, that is) asm,m oo, but as we have seen, in 
return for the decreased freedom, we gain a large number of new properties. 

Elliptic bigrids of a given shape are parametrized by sextuples (a,b;q,t;C,(f>), where C is a genus 1 curve 
over a field k, 4> : C — > P 1 (fc) is a degree 2 function (with associated divisor class r € Pic 2 (C)), q,t <E Pic°(C), 
and a,b e C. The corresponding elliptic bigrid of shape m n is then defined by the formula 

l(0,i,j) = <Kaq i t n - i ) (7.1) 
1 (l,i,j) = cj>(bq-h i - 1 ). (7.2) 

Theorem 7.1. Elliptic bigrids are perfect. 

Proof. Fix a pair (C, cf>), and let the other parameters be generic; also choose a generic point v e C. For each 
partition A C to™, define a BC n (T)-symmetric function f\ by: 

*(•■■*<•■■>= n J^-^Il) R > n) <•••*-" b ^ g - t; c > ^ (7 - 3) 

l<i<n 
l<j<m 

Since fx is BC„ (r)-symmetric, it factors through a symmetric rational function on ¥ 1 (k); consideration of poles 
shows that this rational function is in fact a polynomial. Now, from the vanishing properties of R*£ n \ it follows 
that 

A(...at n -ty"...) = (7.4) 

for \x (f_ A. In particular, for v C to™ with ^ <£ A, we can replace the parts equal to to by any partition with 
parts > to. Thus 

f x (... 1 +(0,i,u i )...) = (7.5) 

for a Zariski-dense collection of extensions of 7. Complementation symmetry of R*^> gives the complementary 
vanishing conditions. In other words, the polynomial associated to f\ is an interpolation polynomial, and thus 
the generic elliptic bigrid is perfect. Since perfection is closed, the theorem follows. □ 
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If two sextuples (a,b;q,t;C,(j)), (a', b'; q' , if; C , <j)') are equivalent in the sense that there exists an isomor- 
phism ip : C — ► C such that 

#»)=«'. <l>(b)=b', #?)=</, f/>{t)=lf, <f>' o ^ = <j>, (7.6) 

the associated elliptic bigrid is the same. In other words, we obtain a map from M. to the space of perfect 
bigrids, where M. is the moduli problem classifying sextuples modulo equivalence. 

Theorem 7.2. Ifm,n > 2, (to, n) ^ (2,2), the map from A4 to the space of bigrids is birational; that is, an 
elliptic bigrid generically determines a unique (up to equivalence) sextuple (a,b;q,t;C,r). Moreover, the closure 
of the space of elliptic bigrids of shape to™ is a rational variety of dimension 8. 

Proof. Let 7 be a generic elliptic bigrid; we need to prove that it determines a unique sextuple. It suffices to 
consider the cases (to, n) € {(3, 2), (2, 3)}, since a larger elliptic bigrid necessarily contains one of the two smaller 
bigrids. Consider the first case. We claim that the sextuple, and thus the full bigrid, is uniquely determined by 
the eight points 

7 (0, 1, 0), 7(0, 1, 1), 7 (0, 1, 2), 7 (0, 2, 0), 7(0, 2, 1), 7(0, 2, 2), 7(0, 1, 1), 7(0, f , 2). (7.7) 

To see this, consider the five pairs of points 

( 7 (0, 1,0), 7 (0, 1,1)), (7(0, 1,1),7(0, 1,2)), ( 7 (0, 2,0), 7 (0, 2,1)), ( 7 (0,2, f ), 7 (0, 2, 2)), ( 7 (1, 1, 1), 7 (1, f , 2)). 

(7.8) 

Each of these is of the form (cf>(x), 4>(qx)) for x 6 C. 

Lemma 7.3. Let C be a genus 1 curve, with (f> : C — > P 1 a degree 2 function, and q G Pic a generic point. 
Then there exists a unique (up to scale) homogeneous polynomial p on P 1 x P 1 of bidegree (2, 2) that vanishes on 
precisely those points of the form {4>(x), 4>{qx)); the resulting subvariety 0/P 1 x P 1 is isomorphic to C . Moreover, 
this polynomial is symmetric under exchanging the factors P 1 . 

Proof. Let <f> have polar divisor (z 1 )(z 2 )- Then evaluation at (<f>(x), 4>(qx)) maps the space of polynomials of 
bidegree (2,2), of dimension 9, to the space C(({zi){z 2 ){zi/ q){z 2 / q)) 2 ), of dimension 8; existence of p follows. 
Any such polynomial cuts out a one-dimensional subscheme of arithmetic genus (2 — 1)(2 — 1) = 1, which is 
therefore the full collection of points (4>(x), 4>(qx)), a curve isomorphic to C. In particular, uniqueness of p 
follows from the connectedness of the subscheme. Since 

{4>(x),<j>(qx)) = (ct>(T/x),<fi(T/qx)), (7.9) 

the subscheme is symmetric, and thus the symmetry of p follows. □ 

Remark. Such "symmetric biquadratic relations" appear in the work of Baxter on the eight-vertex model, and 
other related statistical mechanical models [3]; they have also been mentioned in the context of univariate elliptic 
biorthogonal functions |23|2HI. In particular, the latter reference derives such a relation from the existence of 
an appropriate difference operator. 
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Now, the space of symmetric polynomials of bidegree (2, 2) is six-dimensional, and thus five points in P 1 x P 1 
generically determine a unique (up to scale) such polynomial; a simple computation of minors shows that this 
remains true for our somewhat constrained points. We thus obtain a genus 1 curve Co, together with a degree 
2 map (f>o given by projection onto the first factor. The given five points are then of the form aoio, aogoio, o-o, 
aogo, bo/qQ for uniquely determined points ao, &o € Co, go, to € Pic°(Co). But then 7 is precisely the elliptic 
bigrid corresponding to the sextuple (a , b ; q , t a ; C , <fro), and the theorem follows for (m, n) — (3, 2); the other 
case follows symmetrically. □ 

Remark. In particular, it follows that an elliptic bigrid takes values over a field k if and only if its parameters 
can be defined over that field. 

A similar argument shows that for (m, n) € {(I, 1), (1, 1)} with I > 4, the corresponding map from quintuples 
(forgetting t or q as appropriate) is birational, with rational image of dimension 7; for I < 4, the map is 
surjective, and the range too small for birationality. In the remaining case (m,n) = (2,2), although the space 
of bigrids is 8-dimensional, the map fails to be birational. 

Lemma 7.4. The closure of the space of elliptic bigrids of shape 2 2 is given by the equation 



where \ is the cross-ratio function on P 1 . In particular, this closure is a seven- dimensional rational variety. 
Proof. The four points at, r/ag, 0/q 2 , qr/bt e C multiply to r 2 /g 2 ; it follows that the four points 



(7(0, 1,0), 7(0, 1,1)), (7(0, 2,1), 7(0, 2,0)), (7(1, 1,2), 7(1, 1,1)), (7(1, 2, 1),7(1, 2, 2)) eP 1 x P 1 (7.11) 



satisfy a common bilinear identity (the space of bilinear polynomials and the space C(D) with (D) = r 2 /q 2 
are both 4-dimensional) . The given identity follows. That the image of the space of elliptic bigrids is seven- 
dimensional (and thus agrees with the subscheme cut out by the equation) follows from a direct computation 



While we have so far been unable to derive a full classification result for perfect bigrids-that is, a description 
of its decomposition into irreducible components-we can prove the following partial result. 

Theorem 7.5. Ifm,n> 2, then the space of elliptic bigrids is a component of the space of perfect bigrids; that 
is, the perfect bigrids in a sufficiently small neighborhood of a generic elliptic bigrid are all elliptic. 

Proof. It suffices to show that the tangent space to the space of perfect bigrids at a generic elliptic bigrid has 
the same dimension as the space of elliptic bigrids; that is, that it is seven-dimensional for (m, n) = (2,2), 
and otherwise eight-dimensional. For the small cases (2,2), (2,3), (2,4), (3,2), (4,2), (3,3), we can verify this 
as follows (again a computation facilitated by using the I\ degeneration). First, choose a random instance 
defined over Q, and verify that the tangent space has the required dimension. Moreover, the ideal generated 
by appropriate minors determines a finite set of primes; reducing modulo a prime outside that set verifies the 
result in that characteristic. Then for each remaining characteristic, we can choose a random instance over a 
field of that characteristic, and again verify the dimension. This proves the result in those cases. 



X ( 7 (0, 1,0), 7(0,2, 1), 7 (1, 1,2), 7 (1,2, 1)) = X ( 7 (0,1,1), 7 (0,2,0),7(1, 1,1), 7(1,2,2)) 



(7.10) 



for the 1\ degeneration (see below). 



□ 
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We now proceed by induction; let (m, n) be a pair not on the above list such that the theorem holds for 
all smaller pairs. If m > 4, the two natural sub-bigrids of shape (to — 1)" overlap in a sub-bigrid of shape 
(to — 2) n , and thus by induction, for any bigrid in a small neighborhood of the original bigrid, these sub-bigrids 
are elliptic. By assumption, either to > 5 or n > 3, and thus a generic elliptic bigrid of shape (to — 2)™ uniquely 
determines its parameters. It follows that the sub-bigrids of shape (m — l) n have compatible parameters, and 
thus join to form an elliptic bigrid of shape to" as required. The proof for n > 4 is analogous, and thus the 
theorem follows. □ 

In particular, it follows that, just as we needed to generalize Okounkov's notion of interpolation polynomial 
to obtain elliptic special functions, any attempt to further generalize the elliptic theory must necessarily involve 
a further generalization of balanced interpolation polynomials; there is no room in the space of perfect bigrids 
to add more parameters to elliptic bigrids. 

Now that we have ruled out generalization, it is natural to turn to specialization. The space of elliptic 
bigrids is clearly not closed; it is therefore of interest to describe the various degenerations of elliptic bigrids, 
i.e., to give explicit descriptions of the points in the closure of the space. Again, we will be content with a 
partial result, a description of the four (over the algebraic closure) particular classes of degenerations that can 
produce arbitrarily large regular bigrids. 

The key to understanding these degenerations is the curve Co C P 1 x P 1 constructed in proving Theorem 
17.21 If this curve is smooth, then the bigrid is elliptic, so assume it is singular; we also assume it is reduced 
(since the nonreduced case can only be regular if to = 2, and corresponds to the case q 2 — 1). Thus Co has 
isolated singular points, and since it has arithmetic genus one, consists of a collection of P 1 s. 

Suppose first that Co is irreducible, and thus has only one singular point over k. The singular point must 
lie on the diagonal of P 1 x P 1 , and by simultaneous (fc-rational) linear fractional transformation can be taken 
to be (oo, oo). We thus obtain a curve on A 1 x A 1 = A 2 of the form 



If «20 = 0, any sequence of points xq, x\, X2, X3 € P 1 with (xq, Xx), (xi, X2), {x2, 13) 6 C will contain the same 
point twice; it again follows that this case cannot arise from large regular bigrids. We may thus assume 020 = 1. 
If an 7^ —2, then we can eliminate the linear term; then aoo 7^ 0, lest (0, 0) be another singular point. There 
are then two cases, depending on whether the quadratic term has roots over k. If it does, let q be such that 
q 2 + an g + 1 = 0. Then we find that the smooth fc-rational points on C are precisely those points of the form 



with u € k* . In other words, the smooth points of Co can be identified with the multiplicative group in such 
a way that translation by go is multiplication by q. Similarly, the translation by to induced by the bigrid is 
multiplication by some element t € k* , as is the translation by b/a. We thus obtain the following class of 
(perfect) bigrids, upon restoring the linear fractional transformation freedom. 

Definition 17. A degenerate elliptic bigrid of type 1\ is a bigrid given by the formula 



a 2 o(x 2 + y 2 ) + anxy + a w (x + y) + a 00 . 



(7.12) 



u + a o/u qu + a o I QU 



(7.13) 



q - 1/q q - l/q 



7 (0, = cj>(aqn n ~*) 7 (1, = cfiibq-H 1 - 1 ), 



(7.14) 
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where a, b, q, t e k* , and is a fc-rational degree 2 function taking the same values at and oo. 

Remark 1. In particular, it is quite straightforward to compute the tangent space to a generic such bigrid 
of shape 2 2 , and thus show that that space is 7-dimensional; since the closure of the space of elliptic bigrids 
contains this space, its dimension must be at least as large. 

Remark 2. The notation Ii is motivated by the fact that the curve C is a degenerate elliptic curve of Kodaira 
symbol Ii, and similarly for the labels we use for the other untwisted classes. 

We similarly obtain a twisted version of this, when x 2 + auxy + y 2 is irreducible over k. Let I be the splitting 
field of this polynomial. 

Definition 18. A degenerate elliptic bigrid of type I[ is a bigrid given by the formula 

7(0, - ^{aqH^) 7 (1, = cPibq-i?- 1 ), (7.15) 

where a, b, q, t are in the subgroup of I* of norm 1 and is a nonconstant function on this subgroup of the 
form 

^ = 7 Tr(,) + ( T (? - 16) 
a, (3, 7, 5 e k. Note that 4> takes values in k U {oo}. 

The final case with a single singular point is when the quadratic term has a multiple root; i.e., the singular 
point is a cusp. In finite characteristic p, we find that the translation map has period p (period 4 in characteristic 
2), and thus only leads to bigrids of bounded size; we may thus assume k of characteristic 0. We obtain the 
following. 

Definition 19. A degenerate elliptic bigrid of type II is a bigrid given by the formula 

rf(0,i,j) = <f>(a + jq + (n-i)t) r,(l,i,j) = <t>(b-jq + (i-l)t), (7.17) 

where a, b, q, t £ k, and is a degree two function on k of the form 

, ax(r - x) + (3 ,_ , 

^ = 7 x(r-x) + r (7 ' 18) 

with a, P, 7 , S, r € k. 

Remark. This case can also be thought of as a limit of the I\ case, by taking all of the points to 1 at comparable 
rates, including the zeros and poles of (f>. It corresponds (for sufficiently large to, n) to a 6-dimensional subvariety 
of the space of elliptic grids. 

There are two cases to consider when Co is reducible. Again, to produce regular bigrids, each factor must 
have bidegree at least (1, 1), so there are two factors, both of bidegree (1, 1). Either the individual factors are 
symmetrical, and distinct, or the two factors form an orbit under the exchange map. Note that since the curve 
possesses a large number of rational points, each component must in fact be rational. 

In the second case, we find that either there are two points of intersection (on the diagonal), possibly defined 
over a quadratic extension, or there is one rational point of intersection where the two components are tangent. 
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(I.e., the curve has Kodaira symbol I 2 or III respectively.) For I 2 , we can transform to obtain the curve 
(x — qy){y — qx), while for III we can transform to (x — y) 2 — q 2 . 

Now, suppose the curve associated to j translation is reducible, and consider the curve associated to i 
translation. We find that it must also be reducible, and further that there are compatibility constraints. 
Further analysis reveals the following possibilities for bigrids, covering all of the reducible cases. 

Choose an element k G k* , and construct a group k* .2(k) on the set k* x Z 2 by the rule 

(x,0)(y,0) = (xy,0) (x, 0)(y, 1) = (xy, 1) (x, l)(y,0) = (xy, 1) (x, l)(y, 1) = (xy/ K ,0). (7.19) 
Definition 20. A degenerate elliptic bigrid of type I 2 is a bigrid given by the formula 

7(0,1, j) = 0(aqH n - 1 ) 7 (1, = Mbq-H*- 1 ), (7.20) 
where a, b, q, t € k*.2(n), and <fi is defined by 

0(x,O)=^(r/ a!> l) = ^±| > (7.21) 
7X + 

for some nonconstant, rational, linear fractional transformation and some r G fc*. 

Remark 1. Similarly, there is a twisted version with A;* replaced by the norm 1 subgroup of a quadratic extension, 
and (f> replaced by a nonconstant k U {oo}-valued map 

<K^)=^(r/x,l) = ^±^. (7.22) 

■jx + 7 

Remark 2. In truth, this is comprised of a total of 8 different cases, depending on in which coset of k* each of 
a/6, q, t lies. The result, for sufficiently large m, n, is a 7-dimensional space of bigrids, unless q, £ € fc*, in which 
case the space is 6-dimensional. In particular, only the latter two cases (in which the full group structure is not 
being used) can be obtained by degeneration of I\ . (The corresponding biorthogonal functions are multivariate 
analogues of the biorthogonal functions of pp.) Apparently the other I 2 cases have never been studied, even 
at the univariate level (where for q ^ k* we obtain series such that the term ratios t2ml^2m-\ and t2 m +i/t2m 
are different rational functions of q 2m ). One can thus view such a series as a sum of two basic hypergeometric 
series (by separating out the terms of even index from those of odd index). For instance, one such identity (a 
limiting case of the Frenkel-Turaev summation identity [S]) is that 
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a, bi,b 2 , 63, 64, a 2 /6i6 2 6 3 6 4 
a/61, a/6 2 , a/63, a/6 4 , g 2 6i6 2 6 3 64/a 



2 2 



(a - 6i)(a - 6 2 )(a - 6 3 )(a - 64) 6 5 \q 2 a/bi, q 2 a/b 2 , g 2 a/6 3 , g 2 a/6 4 , g 2 6i6 2 6 3 64/a 



a 2 (l - 6i)(l - 6 2 )(1 - 6 3 )(1 - 64) / a,a 2 6i,g 2 6 2 ,g 2 6 2 ,g 2 63,a 2 64,a 2 /6i6 2 6 3 64 „ , 

q ,q I 

^ / 

(a, a/6i6 2 , a/6163, a/6164, a/6 2 6 3 , a/6 2 6 4 , a/6364, a/6i6 2 6 3 6 4 ; q 2 ) 



(a/61, a/6 2 , a/63, a/64, a/6i6 2 6 3 , a/6i6 2 6 4 , a/616364, a/6 2 6 3 6 4 ; q 2 ) 

as long as both sums terminate. Note that the above hypergeometric sums are balanced but not quite well- 
poised. 
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Remark 3. This is superficially similar to the degeneration of ordinary interpolation polynomials denoted IV 
in |13| . which exists only for n = 2; the latter is not, however, a degeneration of our I2, but rather corresponds 
to a degenerate curve of Kodaira symbol 1^ (which never produces a regular bigrid). 

Finally, we have the additive analogue of type I2] we define k.2(n) analogously. 
Definition 21. A degenerate elliptic bigrid of type III is a bigrid given by the formula 

7(P,i,j)=<Ka+jq+(n-i)t) y(l,i,j)=<t>(b-jq+(i-l)t), (7.23) 
where a, b, q, t S fc.2(/t), and 4> is given by 

^,0)=^(t-i 1 1)=^ (7.24) 
72; + 

for some nonconstant linear fractional transformation. 

Remark. Again, aside from degenerations of type II, and aside from a brief mention of the version of this 
degeneration for ordinary interpolation polynomials in |13| (there denoted by Ill(abc)), this case does not 
appear to have been studied. Each III case is one dimension smaller than the corresponding I2 case. 



8 The trigonometric degeneration 

Of course, the classification of degenerate elliptic bigrids is just a first step towards an understanding of the 
various possible limit cases of our results. The "trigonometric" limit (corresponding to bigrids of type I{) is of 
particular interest, as this is where the Koornwinder polynomials |5] and Okounkov's interpolation polynomials 
(and thus the results of are to be found. Although we will only be considering those degenerations 
necessary to connect our results to those of ]Tjg\, it is striking how many different limits arise there. We thus 
second Rosengren's call in ^7] for a more systematic classification of limits of elliptic hypergeometric structures. 

The only subtle point in extending the algebraic construction to the I\ case is the construction of the 
difference operator; the point is that the divisor (r 2 ) on a smooth C degenerates to a divisor that hits the node. 
As a result, a simple condition on poles is not quite enough to specify 8^; we must add a condition comparing 
the asymptotics along the two branches of the node. We thus obtain the following expression for the degenerate 
( j( n ) j where t n ~ 1 UQUiU2U 3 = r 2 , and all parameters lie in k* . 

S (n) (. . . Xi . . . ;uo, tii, u 2 , u 3 ; t; G m , t) (8.1) 
yr u a t n ~ l (l - u Xt/V)(l - uiXi/r)(l - u 2 Xi/r)(l - u^Xi/r) -j-j- 1-txiXj/r 

Xi(l - X 2 Jt){\ - U 1 l to t"-Vr)(l - U 2 U t n -yT)(l - U 3 U Q t n -*/T) 11 1 - XiXj/T ' 



(Here G m refers to the multiplicative group scheme, or rather its natural compactification to a nodal P 1 .) Then, 
if we define Am (uq; q; G m , r) analogously to the elliptic case, we have the following analogue of Theorem 16.31 

Theorem 8.1. Define a difference operator on BC n (r)- symmetric rational functions by 

(B (n) (uo, ui, u 2 ; q, t;G m , r)/)(. . . x t . . . ) (8.2) 



= IT (1 + t-l T )5 {n) (■ ■ ■ Xi . . . ; tto, ui,u 2 , — s T ;t;G m ,r)f(...x l ...), 

- L - L t n L Unll-\U2 



L U UiU 2 

l<i<n 



■50 



where l t (x) = t/x. Then for any integer m > 0, 

fl (n) ("o, mi, «25 g, t; G m , r)4"' (uo; g; G m , r/q) c 4")("o; g; Gm, r). (8.3) 

Moreover, 

{D (n+k) (u , ui,u 2 ; q, t; G m , r)f)(zi, z n ,u , £ fc_1 u ) (8.4) 

= (L)( n \u , ui,u 2 ;q,t; G m , r)/( , uo, ■ ■ ■ ,i fc_1 «o))(2i, ■ •• ,z n )- 

Proof. The only way in which the proof differs is that we need to show that the image function is smooth at 
the node. This is true by assumption for /, as is the condition (by t r -invariance) that either branch of the 
node gives the same limit. But then multiplying by 5^ gives a function with linear growth at the node, in 
opposite directions along the two branches. Thus symmetrizing kills off this linear growth, giving smoothness 
as required. □ 

Similarly, the A , A and diagonal binomial coefficient symbols are all straightforward to define. We find, if 

U Ul ■ ■ ■ W2TO-1 = t"\ 



U)(x:Uq, ...,U2m-l,G m ,T) = TT 1-, (8.5) 

± 1 1 — U r X T 

0<r<2-m ' 

and thus if wo^i ' ' ' v 2m-i — q m a m , 

Al(a\v ,...,v 2m . 1 ;q,t;G m )= J[ J] jlf"!^? (8-6) 

(i,i)eA0<r<2m rH 1 

Cx(vo,...,v 2m -i;q,t) 



Cl(qa/v , ■ ■ ■ ,qa/v 2m -i;q,t) 

Similarly, if vqv x ■ ■ ■ v 2m -i = {t/q)(qa) m ^ 1 



(8.7) 



A\(a\v , . . . ,v 2m -i;q,t; G m ) = 



q \x\ t 2n(x )c o x2 {qa . q ^ t )C° x (v , . . . , v 2m -i\ q, t) 



C x (q,t;q,t)C^(a,qa/t;q : t)C < Kqa/v ,...,qa/v 2m -i;q,t) 
Finally, 

A\ = bWC+(a;q,t)C° x (l/b,qa/b; q ,t) g) 

[a,b];q,t;G m C \ ( a / b l 1,t)C^(b, qa; q,t) 
(Here we use the C* notations of which are simply the limits p — > of the elliptic notations.) 

It is then trivial to state the analogues of the elliptic results for this most general trigonometric case; 
simply replace C by G m in each identity of Section |SJ This in particular gives multivariate analogues of the 
1009 biorthogonal rational functions of Rahman |14| (or, more precisely, since we are considering only discrete 
biorthogonality here, they are multivariate analogues of the biorthogonal rational functions of Wilson |3U|: in 
any event, continuous biorthogonality is preserved in the trigonometric limit). The new thing that appears is 
that we can take limits as various parameters approach or infinity, i.e., the node of the curve. There are, it 
turns out, three different limits of interpolation functions and binomial coefficients that need to be considered 
for a full understanding of 16 from our current perspective. 

Of course, one of these limits is simply the ordinary interpolation polynomials P*^ (; t, s). More precisely, 
we have the following limit theorem. 
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Theorem 8.2. The interpolation functions and ordinary interpolation polynomials are related by the identity 
p*(n) ( . s t^WcK^^r^av^a/v-g^) (n) 2 

Moreover, replacing the limit b — > oo 6?/ — > gives i/ie same result. 
Proof. Fix a, g, £, and let 

/i n) (...^... ;U ) (8.11) 

denote the corresponding right-hand side. We claim that for each n, (; i>) is a well-defined PC n -symmetric 
Laurent polynomial of degree at most |A|. Indeed, from the branching rule, we find 

^ ^Cl/At n ~ 1 axt l ]<l,t) /A\ , ^ 

„ °A/ K V r « w >?> r J W [t"-ia/6,t]( 9 a/tfc,gt"- 1 D/6.t"- 1 a/D);g.t;G m 

(8.12) 

so by induction it suffices to show that the binomial coefficient has a well-defined limit. This then follows readily 
from the algebraic analogue of Corollary 14. 81 

Thus ji (; i>) is a _BC„-symmetric Laurent polynomial of at most the correct degree, vanishing at points of 
the form . . . at n ~ % q ,li . . . for jj, <f_ A, and is thus a multiple of the ordinary interpolation polynomial as required. 
That this multiple is 1 follows by comparing values at the point . . . vt n ~ l .... 

The argument for b — > is analogous. □ 

Remark 1. One can also use the point . . .at n ~' l q Xi ... to set the normalization; the resulting computation is 
somewhat more complicated, but there is the slight advantage that the proof in ^1 of the latter normalization 
does not use any properties of Macdonald polynomials. As a result, we see that we can deduce such properties 
from the elliptic theory. 

Remark 2. It ought to be possible to prove directly, using only the vanishing conditions, that the interpolation 
polynomials are limits of interpolation functions. More precisely, given any partial order on partitions refining 
the inclusion partial order, the complementary vanishing conditions define a corresponding filtration on the 
space of _BC„-symmetric rational functions; the main question is then to show that one of these filtrations has 
the correct limit as b — > 0,oo. (Two cases of particular interest are the dominance order, for which the limiting 
filtration is simply that given by dominance of monomials, so in particular is independent of q and t, and the 
inclusion order itself, for which the limiting filtration is, by results of ^B], that induced by the Macdonald 
polynomials.) 

Plugging into the respective definitions of binomial coefficients gives the following result (in which the 
right-hand side uses the notation of Section 4 of jTT)]L 



Corollary 8.3. We have the following identity of binomial coefficients. 

(8.13) 



'AN (-l)Mf*00CH-(a;g,i) 
hm 11 



6^0,00 \jmJ [aM . q>t . Gm g»0*'>CJ>(?o; q, t) 



q,t,\/a 



Now, it then follows that we have a corresponding limit for the inverse binomial coefficients of |16| . 
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Corollary 8.4. We have the following identity of binomial coefficients. 

■m ( x ) = Hg!MMW (8 . 14) 

Now, the limiting bigrid on the left is in fact a regular bigrid of type I2', we thus obtain a direct description 
of the inverse binomial coefficients in terms of vanishing conditions. In particular, this explains what happened 
to the evaluation symmetry condition when extending from shifted Macdonald polynomials 20, 8 to ordinary 
interpolation polynomials 12 : rather than relate two functions of the same type, the symmetry relates functions 
from two different degenerations. The resulting family of degenerate interpolation functions is simply the limit 
as r — > or co of the general trigonometric interpolation functions. In particular, for these functions, the 
generalized eigenvalue problem simplifies to an eigenvalue problem, relative to the family of difference operators 
(parametrized by c with a, 6, q, t fixed): 

V t |/|(|/| - 1)/2 IL g /Ml - g/giXl - b/ Xi )) n^ 7 (c(l - Si/cXl - r-iab/cxj)) 
, , 2 r„, IWn(c(l - at n -yc)(l ~ 6t»-*/c)) 

II — — /(• • •'/ 1. (8.15) 

X 7 X •> 

with eigenvalues 

q-WC°(a q t^/c, q c/b;q,t) 

C° x (t^a/c,c/b;q,t;p) ' 1 " > 

(Replace c by r/c in the difference equation for trigonometric interpolation functions, then take r — * 0, 00. Note 
that the further limit c — > 0, 00 makes this the identity operator.) 

If we attempt to obtain the bulk branching rule for ordinary interpolation polynomials (Theorem 3.9 of |16p 
as a limit of our theory, we find a third limit of binomial coefficients appearing. 

Theorem 8.5. We have the following identity of binomial coefficients. 

a^O.oo V/x/ [o ,6]. gi i. Gm t n WC t ,(t; q ,t)C x (b;q,t) /M l-i fc 

Proof. Take the appropriate limit of the bulk branching rule for ordinary interpolation functions to obtain a 
bulk branching rule for ordinary interpolation polynomials. Taking a further limit produces the Macdonald 
polynomials (i.e., polynomials satisfying Macdonald's difference equation); the result then follows essentially by 
the definition of P\//i- ^ 

Remark 1. If we take a further limit 6 — > 0, the binomial coefficient can be expressed in shifted Macdonald 
polynomials; we thus obtain Theorem 1.4 of [7], which expresses the corresponding binomial coefficients via the 
plethystic expression 

Pv^t^]; <?,*)■ (8.18) 

Remark 2. Of our three types of degenerate trigonometric binomial coefficients, this is the only one that includes 
the four important special cases 

'AN AN /AN /A N 



(8.19) 

/V [a, !]:<?, t \Av [a,l/q];q,t \Av [a,t];q,t \^/ [o,g/t];g,i 
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That the first three specializations of skew Macdonald polynomials have nice factorizations is well known (the 
first from the fact that P\/^ has positive degree, so constant term 0, the second and third from the coefficients 
of the Pieri identities and branching rule). That 



1 - (a/t) k 

p a/„Q -tiy ];<?,*) (8- 2 o) 



can be expressed as a product of binomials does not appear to have been observed before. 

These also correspond to a degenerate elliptic bigrid of type I2, this time with a/b ^ k* . The difference 
operators are the same as for the previous I2 case, except that now the family has a, c, q, t fixed and b variable; 
in particular, we return to a generalized eigenvalue problem scenario. 

We thus find that all of the ordinary interpolation polynomial identities of are indeed limits of identities 
satisfied by our interpolation functions. Similarly, the Koornwinder polynomials are limiting cases of our 
biorthogonal functions: 

(")/ ^ ■ + ■+ + + „/4.2n-2. . n /i\2\_ K \ (■ ■ ■ x i ■ ■ ■ i 3, *! *0, *1, *2, £3) 



lim R [n) (. . .Xi . . . ;t :ti,t 2 ,t 3 ;uo,q/t " £0*1*2*3^0; Q, *; G m , (1) ) = 



u °^°° (. . . tot"-' . . . ; q, t; t , ti, t 2 , t 3 ) 

(8.21) 

where the normalization factor on the right can be deduced by taking coefficients of PT^ n ' (;q,t, to) on both 
sides. In particular, we note that Theorems 15.61 15.91 and 15.101 give rise to new results at the Koornwinder 
level, generalizing similar results from |16j in which the relevant trigonometric binomial coefficient has second 
parameter E {l/q,t} (i.e., corresponds to a difference or integral operator). 

As an indication of some of the subtle issues that can arise when degenerating the elliptic theory, we 
consider the following "rational" limit. Here we take q, t — * 1 along a one-parameter family in such a way that 
the corresponding maps of tangent spaces have ratio a at the limit point. In other words, we blow up Pic°(C) 2 
at the point (1, 1) and consider the interpolation function at a point on the exceptional divisor. Note that the 
result is essentially independent of every parameter except a, including the choice of C itself. 

Theorem 8.6. Let (C, t) be an arbitrary genus 1 hyperelliptic curve over a field of characteristic 0. Then we 
have the following identity of interpolation functions, for a, b, v S C , a £ k* generic. 

r D *(n)/ , / s . n s J\{. ■ -x{zi,a,v,b\C, r) . . . ; a) 

hm R x y ' {...z t ...;a,b (v );q,t;C,T) = — , 8.22) 

j,t->i J A (. ..l...;a) 

aq/at — >a 

where Jx(](x) is a Jack polynomial. 

Proof. It suffices to prove the result for (C, r) of the form (C*/p, (l) 2 ), in which case taking a limit in the 
branching rule gives the desired result. □ 

There are several points to observe here. The first is that, unlike in the case of bigrids, we obtained a 
nontrivial degeneration without degenerating the curve itself. In fact, while the corresponding elliptic bigrid has 
a well-defined limit, that limit is very far from being regular; indeed, the limiting bigrid is actually independent 
of i and j, having two constant values independent of a. In other words, the map from elliptic bigrids to 
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interpolation polynomials is merely birational; a full classification of degenerate cases of the one is insufficient 
to give a classification the other. The second point is that, compared to the Macdonald polynomial case, it 
is much harder to obtain a useful limit from the difference equation. Indeed, the difference operator actually 
becomes the identity in the limit, and even the generalized eigenvalue problem, while it has a nontrivial limit 
when properly scaled, simply gives (after a nontrivial calculation) the equation 



corresponding to the fact that J\ is homogeneous of degree |A|. The third point is that we can similarly obtain 
limits of binomial coefficients analogous to the three trigonometric limits we obtained above; we find again that 
those limits are essentially independent of those parameters not approaching 1. We finally note that using the 
binomial formula, we can obtain an analogous limit for the biorthogonal functions; if toti, £2^3, uo u i all tend to 
I simultaneously with q, t, then the biorthogonal function becomes a multivariate Jacobi polynomial, up to an 
analogous change of variables. In the case to ~ I, t% ~ —I, the continuous biorthogonality density of ^Hl easily 
simplifies to the usual Jacobi ensemble (just as with the analogous limit for the Koornwinder polynomials); 
more generally, if \to\ = 1 = | ^2 1 , we can at least formally obtain the appropriate ensemble, using the remark 
following Lemma 3.3 of |15|. 

9 Open problems 

It seems fitting to conclude by discussing several open problems arising from the above theory. We have 
already discussed partial results for one such problem in the previous section, namely the question of classifying 
degenerations of interpolation and biorthogonal functions; along these lines, we also note the question of finding 
a direct proof that the interpolation polynomials are limits of the interpolation functions. For the remaining 
problems we consider, we have even less in the way of explicit results or even conjectures. 

The most popular approach to the theory of Koornwinder polynomials involves the so-called "double affine 
Hecke algebra" along with the associated representation of the affine Hecke algebra. A major open question, 
therefore, is how to generalize this theory to the elliptic level, both to understand how the Koornwinder theory 
relates to the elliptic theory, and ideally to allow the elliptic theory to be extended to other root systems. 
Via evaluation symmetry, this should be related to the question of finding a general Pieri identity for the 
biorthogonal functions. We note the following partial result. 

Proposition 9.1. For otherwise generic parameters satisfying t tot\t2t^UoUx = qr 3 , and generic C , 

let c\ K be the coefficient of 




(8.23) 



R^ l \. ..Zi...\ t :t\,h,h\ u q, Ut/q; q, t; C, r) 



(9.1) 



in the product 



{ I x(zi,u,vt n \u ]C,T)R^\...Zi...]t :ti,t 2 ,h]u ,ui]C,T). 



(9.2) 



Ki<n 



Then c\ K = unless A C 1" + n C 2™ + A. 
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Proof. That k C 1™ + A follows from the Pieri identity for interpolation functions and the fact that the 
biorthogonal functions are triangular in the appropriate basis of interpolation functions. That Ac 1" + k 
follows by the fact that for t" -1 q m to£i = t, c\ k can be computed via biorthogonality, and thus satisfies a 
symmetry switching A and n. □ 

In the cases u £ {to,ti,t2,t3,to/q,t\/q,t2/q,t3/q}, we can combine the quasi-Pieri identity, Theorem 15. 101 
with the connection coefficient identity, to express the above coefficient as a sum over binomial coefficients with 
b = 1/q, which is thus explicitly supported on the stated set of partitions. The question is thus whether a 
corresponding formula exists for general u. (An expression in a double sum is straightforward, but it is then 
nonobvious that the coefficients have the support they do.) This would give a family difference operators via 
evaluation symmetry, which (by consideration of the Koornwinder case) would give the analogue of the degree 
1 subspace of center of the affinc Hccke algebra. (See also the remark following Theorem 8.9 of |15|.) 

Similar considerations hold for a general branching rule. 

Proposition 9.2. For otherwise generic parameters satisfying t 2n ~ 2 totit2t3UQUi = qr 3 , and generic u 6 C , let 
c\ K be the coefficient of 

^ n_1) (- t :tx,t 2 ,t 3 ; u , Ul t 2 ;q, t; C, r) (9.3) 

in the specialization 

f6™\. ..z%... ,u;to:ti,t2,tr,uo,uv,C,T). (9.4) 

Then c\ K = unless k <' 2 X. 

Proof. Using the Cauchy identity, Theorem l5.11l one can express the coefficients of the branching rule in terms 
of the coefficients of the Pieri identity (and vice versa); the stated vanishing property follows. □ 

In particular, we see that it is hopeless to expect anything better than a sum for these coefficients, as even 
in the case n = 1 we obtain a nontrivial hypergeometric sum. We also note that such an expression for the 
branching rule would be a new result even for the Koornwinder case. 

Another question, which is less important, but likely to have interesting consequences, concerns the classi- 
fication of regular perfect bigrids. The four non-elliptic components we have identified above all correspond to 
important special cases of the interpolation functions, and thus of the corresponding hypergeometric identities. 
It thus seems likely that other components would give similarly interesting special cases. One such component 
would appear to correspond to the following. 

Conjecture 1. Any bigrid such that = 7(0, i + whenever both sides are defined is perfect. 

Note that this would give a ((n + l)(m + 1) — 2)-dimensional space of perfect grids, and would thus be 
larger even than the m(n + l)-dimensional Cauchy case. When m = n = 2 (in which case the conjecture is 
straightforward to verify), we obtain a 7-dimensional space, distinct from the 7-dimensional space of elliptic 
bigrids; in contrast, a bigrid of shape 2 2 coming from any of the four known cases is necessarily elliptic. In the 
elliptic case, this corresponds to binomial coefficients with b = 1/t, and thus should correspond in some sense 
to an inverse of the integral operator of ^Hl • 
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The other conjectured component we have found (by computing tangent spaces to random elliptic bigrids 
over small finite fields) is the following. 

Conjecture 2. Let t be a linear fractional transformation of order 2. Then any bigrid satisfying 

7 (0, = r( 7 (l,i,i + 1)) = r( 7 (0,i + 1, j + 1)) (9.5) 
whenever both sides are defined is perfect. 

Finally, we observe that we have only extended about half of the results of ^B] to the elliptic case. The 
remaining results fall into two main classes. The first of these is the construction of families of symmetric 
functions algebraically continuing interpolation and Koornwinder polynomials in the number of variables (via 
the quantity t n ). Although it seems unlikely that these results extend cleanly to the elliptic level (given the 
absence of an appropriate theory of symmetric functions) , there should still be some weak version; after all, the 
connection coefficients that arise above depend on the number of variables only via t n , and similarly for the 
bulk branching rule. The other class of results concerns the "vanishing integrals" conjectured in [T^], which 
generalize quadratic transformations of univariate hypergeometric series; elliptic analogues of these would likely 
be interesting even at the univariate level. Note in particular that there are two potential kinds of quadratic 
transformation at the elliptic level, depending on whether the relevant isogeny is a multiplication map (squaring, 
in our notation), or simply an isogeny of degree 2. 
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